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ABSTRACT 



There are substantial quantum fluctuations even in a pure coherent state. For example, the 
quadrature fluctuation in coherent state |q; > is | and the fluctuation in photon number is jap. 
Further, randomization in the case of general classical states which are random superposition of 
coherent states can only increase these fluctuations. But there are some phenomena in which 
quantum fluctuations reduce below the coherent state level. For example, an electromagnetic 
field is said to be electrically squeezed field if uncertainties in the quadrature phase observable X 
is less than the coherent state level (i.e. ( AX)^ < |). Correspondingly a magnetically squeezed 
field is one for which [AXy < |. On the other hand, antibunching is a phenomenon in which 
the fluctuation in photon number is reduced below the Poisson level (i.e (AiV)^ < < iV >). 
All phenomena in which quantum fluctuations reduce below those of the coherent state level are 
called nonclassical. It can be noted that the interaction of dielectric medium with an intense 
electromagnetic field can gives rise to these nonclassical nonlinear optical effects. These effects 
which help us to understand the nature of quantum world in further detail are also important from 
the application point of view. Keeping that in mind, the present thesis is devoted to the study of 
quantum fluctuations of coherent light in non-linear media. 

The interaction of coherent light with a nonlinear medium is modeled here by a general quan- 
tum anharmonic oscillator. The model of a quartic anharmonic oscillator emerges if the lowest 
order of nonlinearity (i.e third order nonlinear susceptibility) is assumed. The quartic oscilla- 
tor model has considerable importance in the study of non-linear and quantum optical effects 
present in a nonlinear medium of inversion symmetry. Silica crystals constitute an inversion 
symmetric third order nonlinear medium and these crystals are used to construct optical fibers. 
In optical communications electromagnetic beams pass through these fibers and the interaction 
of the electromagnetic field (single mode) with the fiber can be described by the one dimensional 
quartic anharmonic oscillator Hamiltonian. Depending on the nature of nonlinearity in physical 
problems the treatment of higher anharmonic oscillators assumes significance. 

Anharmonic oscillator models are not exactly solvable in a closed analytical form. But we 
need operator solutions of the equations of motion corresponding to these models in order to 
study the quantum fluctuations of coherent light in nonlinear media. So we have two alterna- 
tives, either we can use an approximate Hamiltonian which is exactly solvable or we can use 
an approximate operator solution. Operator solutions (solutions in the Heisenberg approach) of 
anharmonic oscillator problems were not available since the existing methods tend to introduce 
inordinate mathematical complications in a detailed study. Due to the unavailability of the op- 
erator solutions people were bound to use rotating wave approximated Hamiltonians to study 
the quantum fluctuations of coherent light in nonlinear media. The situation has improved in the 
recent past and many solutions of anharmonic oscillators have appeared. Some of these solutions 
are obtained as a part of the present work. For example, we have constructed a second order oper- 
ator solution for quartic oscillator and have generalized all first order operator solutions available 
for the quartic oscillator to the m-th anharmonic oscillator. From the generalized solutions we 
observe that there exists an apparent discrepancy between the solutions obtained by different 



4 



techniques. Then the question arises: Which solution should be used for physical applications. 

Therefore, we compare different solutions and conclude that all correct solutions are equivalent 
and the apparent discrepancy is due to the use of different ordering of the operators. 

We use these solutions to investigate the possibilities of observing different optical phenom- 
ena in a nonlinear dielectric medium. To be precise, we have studied quantum phase fluctuations 
of coherent light in third order inversion symmetric nonlinear medium. Fluctuations in phase 
space quadrature for the same system are studied and the possibility of generating squeezed state 
is reported. Fluctuations in photon number are studied and the nonclassical phenomenon of anti- 
bunching is predicted. We have generalized the results obtained for third order nonlinear media 
and have studied the interaction of an intense laser beam with a general (m — l)-th order nonlin- 
ear medium. Aharonov Anandan nonadiabatic geometric phase is also discussed in the context 
of (m — l)-th order nonlinear medium. 
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Chapter 1 

Interaction of electromagnetic field with 
matter 

Since the invention of laser people have observed different nonlinear optical effects in dielectric 
media interacting with an intense electromagnetic field. Studies into these effects are extremely 
important because some of them (such as squeezing and antibunching of photons) do not have 
any classical analogue. These purely quantum mechanical effects are also important from the 
application point of view. Keeping that in mind, in the present work we have discussed a class 
of physical systems in which different nonlinear and quantum optical effects may be observed. 
To be more specific, we consider a physical system in which an intense electromagnetic field 
is interacting with an (m — l)-th order nonlinear medium. The physical system is modeled by 
a generalized anharmonic oscillator (AHO). The operator solution of the equation of motion 
corresponding to our model Hamiltonian is developed using different techniques and is exploited 
to study the possibilities of observing higher harmonic generation, bunching and anti-bunching 
of photons, self induced transparency and squeezing. In addition to this, fluctuations in quantum 
phase of the output field is also studied with the knowledge of field operators. Geometric phase 
of the output field is also taken care of. 

1.1 The model 

The present thesis is devoted to the study of the interaction of a single-mode intense electro- 
magnetic field with a nonlinear medium. In such a nonlinear interaction nonclassical effects like 
squeezing and antibunching of photons may be produced. In fact, here we study the possibili- 
ties of observing these nonclassical phenomena in the physical system of our interest. We can 
continue our study by using either a semiclassical theory or a quantum theory. Actually, a semi- 
classical theory, where the radiation is represented by a classical wave and the atom is quantized, 
can treat many problems in light-matter interaction. However, there are cases where a semi- 
classical theory turns out to be inadequate. These include spontaneous emission. Lamb shift, 
resonance fluorescence, the anomalous gyromagnetic moment of the electron and "nonclassical" 
states of light (e.g. squeezed states). So if we need to study the possibilities of observing non- 
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classical states then we have to use a quantum theory instead of a semiclassical one and as an 
essential requirement of the quantum theory we have to quantize the electromagnetic field. In the 
next subsection the electromagnetic field is quantized and the model Hamiltonian of our physical 
system is constructed in the subsequent subsection. 



1.1.1 Quantization of the electromagnetic field 

To quantize electromagnetic field let us recall some basic results of electrodynamics. The 
Maxwell equations describing the propagation of electromagnetic fields are 

VxE = -f 
V.D = p 

VxH = J + f ^^-^^ 
V.B = 

where J is the current density; the other symbols have their conventional definitions. In a homo- 
geneous isotropic medium B and D are related to H and E by 

B = /xH 

D = eE ^ ■ ' 

where /i and e are, respectively, the magnetic permeability and the dielectric constant of the 
medium. If we restrict ourselves to a charge free medium, then we have 

P = ^ (13) 

Consider the electric field E and magnetic field H inside a volume V bounded by a surface S 
of perfect conductivity. The tangential component ofE,— nxnxE, and the normal component 
of H, n.H will both be zero on S (n is the unit vector normal to S). We can expand E and H 
in terms of two orthogonal sets of vector fields E^ and Hq, respectively. These sets which were 
originally introduced by Slater ^ and nicely explained by Yariv [0] obey the relations 

ka^a = V X Ha 

and 

kjia = V X E„ (1.5) 

where ka is to be considered, for the moment, a constant. The tangential component of E^ on S 
is zero, i.e 

nxEa = Oori^. (1.6) 
If we take the curl of both sides of the equations ( ]1.4| - p3| ) and use the identity 

V X V X A = V(V.A) - V^A (1.7) 
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then we obtain the familiar wave equations 

V^E, = klEa 



Now we can write the total resonator fields of E and H as 

E(r,t) = -E,^g.(t)E.(r) 

H(r,t) = Ea^.Pam.ir) ^'-^^ 
where cOa = Qait) andpa(t) are measures of the field amplitude in a-th mode. Substituting 



( |0[ ) in the first two Maxwell equations ( [TTT| ) and using ([L^-[T3|) we obtain 



Pa = -^ala 
Qa = Pa- 



From ( Jl.lOt ) we have 

qa + iolqa = 

pa + UJ^Pa = . 

This identifies cUa = -7^ as the radian oscillation frequency of the a-th mode. 



(1.10) 



(1.11) 



/Me 

From ( |L9| ) we can easily show that 



and 



qa{t) = / E{r,t).E^{r)dv (1.12) 

Pait) = ^ |^H(r,t).H,(r)rft;. (1.13) 

Now we conclude that the electromagnetic field can be specified either by E(r, t) and H(r, t) or, 
alternatively, by the dynamical variables qa{t) and Pa{t)- The total energy Hamiltonian is 

Ho = - f (fiH.H + eE.E)dv (1.14) 
2 Jv 

where the subscript zero is used to distinguish the free field Hamiltonian (Hq) from the total 
Hamiltonian of the system. Substituting the expansions of E and H ( |1.9D in ( |1.14] ) we obtain 

^0 = E^(Pa + ^aga). (1-15) 
a ^ 

This has the basic form of a sum of harmonic oscillator Hamiltonians. The dynamical variables 
Pa and qa constitute canonically conjugate variables. This can be seen by considering Hamilton's 
equations of motion relating pa to qa and qa to pa- The equations are 

I I "'''^'^ (1.16) 

la Qp^ Pa ■ 
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These are foraially identical with ID , obtained from Maxwell's equations. 

The quantization of the electromagnetic field is finally achieved by considering pa and 
as formally equivalent to the momentum and coordinate of a quantum mechanical harmonic 
oscillator, thus taking the commutator relations connecting the dynamical variables as 



[ga, qb] = 







where we have chosen to work in a system of units in which fi 
bo sonic creation and annihilation operators as 



(1.17) 

1 . At this point we can introduce 



a\{t) 
ai(t) 



Solving dl.lSD for pi and qi we obtain 

Piit) -- 



(2^) 



2 

(2^) 



[uJiqi{t) -ipi{t)] 
[uJiqi{t) + ipi{t)] . 



a|(t) - ai(t) 
al{t) + ai{t) 



(1.18) 



(1.19) 



Now we can express the free field Hamiltonian in terms of annihilation and creation operators as 



Ho = "^uiiajai + 



(1.20) 



Thus we have quantized the electromagnetic field in vacuum. But the aim of the present thesis is 
to study the interaction of a single-mode intense electromagnetic field with a nonlinear medium. 
Therefore, the interaction part of the Hamiltonian should be taken into account. 



1.1.1.1 Conventions used 



Before we go into further detail let us state the conventions which are used in the present thesis. 

i) We work in units in which h = 1. 

ii) The applied electromagnetic field is chosen to be a single mode electromagnetic field 
having unit frequency (i.e u = 1). 

iii) Mass of the harmonic oscillator equivalent to the applied field is taken to be unity. 

iv) Lower case x{t) and x{t) are to be treated as time development of the position and mo- 
mentum variables in the classical sense and the equivalent operator representation are to be made 
by using the corresponding upper case letters X{t) and X{t). Thus 

1 



q{t) = X{t) 



V2 



(1.21) 



(since m = 1) and 



p(t) = mq{t) = mX{t) = X{t) 



[Xit),Xit)]=i. 



V2 



a' — a 



(1.22) 
(1.23) 
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1.1.2 Construction of the model Hamiltonian 



With the above conventions the free field Hamiltonian ( |1.20| ) of our system reduces to 



^0 - 2 + 2 (1.24) 

a' a + 



2- 



Again from ( |1.9| ) we observe that 



and 



E oc g oc (a"^ + a) (1.25) 

'B (X p oc (^a^ — . (1-26) 

Now an intense electromagnetic field interacting with a dielectric medium induces a macro- 
scopic polarization (P) having a general form 

P = XiE + X2EE + X3EEE + .... (1.27) 

where E is the electric field and xi is the linear susceptibility. The parameters X2 and xs 
second and third order nonlinear susceptibilities respectively. Presence of the dielectric medium 
contributes to the electromagnetic energy density. This contribution is proportional to P.E. 
Therefore, the interaction Hamiltonian Hemi is 

i/e™ = A.P.E (1.28) 

where Ac is the proportionality constant. Now, if the symmetry of the medium is chosen in 
such a way that the only nonlinear interaction in the medium appears due to the presence of the 
(m — l)-th order susceptibility (Xm-i) and the macroscopic magnetization (if any) is neglected, 
then the leading contribution to the interaction part of the electromagnetic energy comes through 
the coupling of the (m — l)-th order nonlinear susceptibility. Therefore, the interaction energy is 
proportional to the m-th power of the electric field. Now using (|1.24| - [L28b the total Hamiltonian 



of the system which is a sum of free field and interaction Hamiltonians, can be written as 



H = a)a + \ -\ ^-m-fa"'" + a 



m 



The parameter A is the coupling constant and is a function of the (m — l)-th order nonlinear 
susceptibility (Xm-i) of the medium. The above Hamiltonian ( |1.29| ) represents a generalized 



anharmonic oscillator of unit mass and unit frequency. The equation of motion corresponding to 



(I1.29D is 

X + X + AX'"-^ = (1.30) 

which can not be solved exactly for m > 2. Present work provides a first order operator so- 
lution for this generalized anharmonic oscillator and the solution is used here to study various 
nonclassical properties of radiation field interacting with a nonlinear medium. 
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1.1.2.1 Special cases of the model Hamiltonian 

The Hamiltonian ( 11.29|) represents a class of physical systems. Among those systems m = 4 
is a case of particular interest because even order susceptibilities {x2, X4 ^tc.) would vanish for 
an inversion symmetric medium. Hence the leading contribution to the nonlinear polarization 
in an inversion symmetric medium comes through the third order susceptibility (xs) ^iid the 
corresponding Hamiltonian is 

2 ^ 2 ^ 4^ (1.31) 

= a^a + j^{a^ + a)^ . 
This is the Hamiltonian of a quartic anharmonic oscillator whose equation of motion is 

X + X + XX^ = 0. (1.32) 



Equation ( |1.32[ ) is very important in the context of quantum optics, nonlinear optics, molecular 



physics, (p'^ field theory and many other branches of physics. But equation ( |1.32D is not ex 



actly solvable due to the presence of the off-diagonal terms in the Hamiltonian ( |1.31| ). These 
off-diagonal interaction terms are called photon number non conserving termsQ. The interac- 
tion involving the operator (a^^) may be viewed as the loss (gain) of two photons. On the 
other hand, interaction involving a^a keeps the photon number conserved (i.e one photon is 
created and one photon is destroyed). So the terms proportional to a^, a)"^ , a) a?, a)'^a, and 
a^'^ are photon number nonconserving in nature. It is usually assumed that the photon number 
nonconserving terms have no significant contribution in the time development of the field op- 
erators. This assumption may be realized in the following way. In the absence of interaction, 
the time development of the annihilation operator a(t) = a(0)e^**. Thus the terms proportional 
to a^, a^^, a^a^, a'^^a, and a)^ are rapidly oscillating compared to the terms a)a and a)'^a?. 
These rapidly oscillating terms contribute little to the interaction Hamiltonian. The assumption is 
widely used and is called the Rotating Wave Approximation (RWA). Under RWA ( |1.31[ ) assumes 
an extremely simple and exactly solvable form 

H = a' a + Xga^ a^. (1.33) 
Gerry and Lynch used this rotating wave approximated Hamiltonian (|1.33| ) to study 



the quantum phase fluctuations of coherent light interacting with a third order nonlinear medium. 
Joshi et al [0] used the Hamiltonian ( |1.33[ ) to study the sensitivity of nonadiabatic geometric 
phase on the initial photon statistics in a dispersive fiber. Tanas M\ used the same Hamiltonian 



( |1.33[ ) to study the possibilities of generating squeezed state in the interaction of single mode 
electromagnetic field with a nonabsorbing nonlinear medium. 

The RWA is hardly valid if the noncatalytic nonlinearities are present in the system [Q. For 
example, the term proportional to a^'^ is responsible for higher harmonic generation and occurs 



' According to Loudon the off-diagonal terms are called energy nonconserving but Gerry [Q] mentioned 
them as photon number nonconserving terms. The phrase 'photon number nonconserving' appears to be more 
appropriate in the context of the physical process involved. So we mention the offdiagonal terms as photon number 
nonconserving terms. 
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naturally corresponding to the field interacting with a nonlinear medium. In fact, the terms 
proportional to a'^ and a^'^ were taken into consideration by Tombesi and Mecozzi [ [T0| ] to study 
the possibilities of squeezing (of coherent light passing through a nonlinear medium). Again the 
nonconserving energy terms are responsible for the well known Bloch-Siegert shift if the field 
frequency is low (e.g RF and MW) [pT]]. 

Keeping these facts in mind, in the present work, we include the offdiagonal terms in the 
model Hamiltonian and observe that important physical information is lost under RWA calcu- 
lations. In spite of working with ^3 medium only, here we study (m — l)-th order nonlinear 
medium in general with particular emphasis to the third order nonlinear medium, i.e. to the case 
of a quartic anharmonic oscillator. 

Depending upon the order of the anharmonicity the anharmonic oscillators are identified by 
different names such as quartic and sextic anharmonic oscillators. Hamiltonian of these oscilla- 
tors, their equations of motion and classical counter parts of the equations of motion are required 
by us at different stages of the development of the present work. So here we give a brief table, 
containing this information. 



Anharmonic 
oscillator 


Classical 
Hamiltonian 


Classical equation 
of motion 


Quantum 
Hamiltonian 


Quantum equation 
of motion 


Nonvanishing 
susceptibility 


Generalized 


2 -2 \ 

2 2 


X + X + Ax™-1 = 




X + X + AX™-! =0 


Xm-l 


Quartic 


if. ^ if. _|_ A-r* 
2 ^ 2 ^ 4-^ 


X + X + Xx^ = 


X^ . X^ . X y4, 


X + X + AX3 = 


X3 


Sextic 


x^ , x^ , \ 6 
2 ^ 2 ^ 6 


X + X + \x^ = 


2^2^ 6 


X + X + AX^ = 


X5 


Octic 


x^ , x'^ , X S 
2 -t- 2 8 


X + X + Xx^ = 


+ + 


X + X + AX^ = 


XT 



Table 1.1 

The model Hamiltonian ( |1.29| ) is not exactly solvable in a closed analytical form. So we have 
two alternatives, either we have to approximate the Hamiltonian itself by neglecting some terms 
(as it is done in RWA) or we have to find approximate solution of the equation of motion ( |1.30D 
corresponding to the Hamiltonian ( [1.291 ). The first option is already tested for xs medium by 
various authors [|]-|]]. Therefore, here we approach the problem by using the second option. The 
approximate solution can be obtained by different approaches. In the following we discuss the 
different possible approaches to the AHO model. 



1.2 Different approaches to the model 

In the last section interaction of an intense electromagnetic field with an (m — l)-th order non- 
linear medium is modeled by a generalized anharmonic oscillator. The equation of motion cor- 
responding to this model Hamiltonian is not exactly solvable and, in-fact, this provides one of 
the simplest examples of quantum mechanical systems which can not be solved without making 
use of approximation techniques. Approximate solution of the AHO problem can be obtained in 
three different pictures in quantum mechanics. These approaches are called Schrodinger picture, 
Heisenberg picture and interaction picture. At first we discuss these approaches briefly and after 
that we choose the appropriate one for the further development of the present study. 
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1.2.1 Schrodinger picture 



The Schrodinger picture deals with the time development of the wave function. Basically, the 
Schrodinger picture is used to solve the quantum oscillators as eigenvalue problems. In these 
problems, energy eigenvalues are expressed as the sum of different orders of anharmonic constant 
A. However, the energy eigenvalues are found to diverge for large anharmonic constant. In case 
of small A, the eigenvalues for different orders can be summed up and the convergence of these 
sums are ensured by the Borel summability [|T2t [T3ll and/or Stieljet conditions [ [l4|l . Schrodinger 



approach is found very successful to treat the quantum anharmonic oscillator problem and there 
exists an extensive amount of literature on this approach [ [T^ - |2^ ]. Schrodinger approach is also 
successfully used to enrich the subject of large order perturbation theory [[T51-[Ty|], divergent ex- 



pansion of quantum mechanics [ [T2| , ^and Laplace transformation representation of energy 
eigenvalues [|TJ] and computational physics [ p2| - p5| ]. 

1.2.2 Heisenberg picture 

In this approach, one solves Heisenberg equation of motion to obtain the time evolution of the 
operators. But the noncommuting nature of the operators poses serious difficulties to obtain 
approximate operator solutions of the quantum anharmonic oscillators. Only few methods are 
available till date. In 1967 Aks [ ]2^ ] obtained an operator solution of a quantum quartic anhar- 
monic oscillator (QQAHO) by using the renormalization technique. In 1970 Aks and Caharat 
[ p9| ] reported an extension of the earlier work by using the method of Bogoliubov and Krylov. 
But the recent interest in the problem started with the work of Bender and Bettencourt [ ]30| - pTl ] 
who obtained the solution of a QQAHO using multiscale perturbation theory (MSPT) in 1996. 



Later on Mandal [ |32| ] proposed a Taylor series approach, Egusquiza and Basagoiti [ [33D devel 



oped renormalization technique, Kahn and Zarmi p^] used a near identity transform method to 



solve the QQAHO problem. Fernandez [ |35|J360 used an eigenvalue approach to solve QQAHO. 

Some of the recently derived solutions are obtained as a part of the present work [^7|-^ and 
are used here to study quantum optical and nonlinear optical effects. Therefore, these solutions 
are much important in the context of the present thesis. Keeping that in mind, here we give a list 
of the solutions of AHO problem available in Heisenberg picture. 
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A lit nofc 




iVlClllvJU 




aks yzpy 


1 OAT 


Renormalization 


Quartic 


Ai^o /^oVioT«o<^ noon 
Aks ana L-anarat n£y\\ 




Bogoliubov and Krylov 


Quartic 


Bender and Bettencourt []3^] 




A /TCDT' 

Mbrl 


Quartic 


Mandal [pzy 


lyy© 


laylor series 


Quartic 


Egusquiza and Basagoiti [33] 


lyyo 


Renormalization Group 


Quartic 


Kahn and Zarmi [|34j] 


lyyy 


Near identity transform 


Quartic 


Speliotopoulos [|4T|] 


OAAA 
ZUUU 


nonperturbative 


Quartic 


Pathak [g7|] 


OAAA 

zUUU 




Generalized 


Pathak and Mandal [|8p 


OAA1 


laylor series 


Quartic 


Pathak and Mandal [Q 




Taylor Series 


Sextic and Octic 


Fernandez [ 




2001 


eigenvalue approach 


Quartic 


Fernandez [ 




2001 


eigenvalue approach 


Particular values of m 


Pathak and Fernandez [|40|] 


2001 


All the existing methods 


Generalized 



Table 1.2 



1.2.3 Interaction picture 

In this approach the time dependence of the physically measurable quantity is described by letting 
both observable and state vary with time. The observables are rotated in one direction with a 
transformation generated by part of the Hamiltonian, and the states are rotated in the opposite 
direction by a transformation generated by the other part (the interaction part) of the Hamiltonian. 
In case of anharmonic oscillator there is a simple relation between the annihilation operator 
obtained in Heisenberg approach and that obtained in interaction picture. Only one solution [ pO| ] 
of general anharmonic oscillator in this picture is reported till date. 



1.3 Quantum fluctuations and appropriate approach for our 
study 

Second order variance {AAY =< > — < A is a measure of quantum fluctuations 
associated with an arbitrary quantum mechanical observable A. In the present work quantum 
fluctuations (variance) are calculated with respect to initial the coherent state \a >. There are 
substantial quantum fluctuations even in a pure coherent state. For example, the quadrature 
fluctuation in |a > is 

(AX)2 = 1 and (AX)2 = 1 (1.34) 
while the fluctuation in photon number is 

(AN)^ =< N >= lal"^ (1.35) 

where = d^a, is the number operator. Further randomization in the case of general classical 
states which are random superposition of coherent states can only increase theses fluctuations. 
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But there are some phenomena in which quantum fluctuations reduce below the coherent state 
level. Let us give some examples, 

Example 1: An electromagnetic field is said to be electrically squeezed field if uncertainties 
in the quadrature phase observable X is less than the coherent state level (i.e. (AX)^ < |). 
Correspondingly a magnetically squeezed field is one for which (AX)^ < |. 
Example 2: In the phenomenon of antibunching (single mode) the fluctuation in photon number 
is reduced below the Poisson level ( |1.35| ), 



[ANf <<N> 



(1.36) 



Example 3: The usual parameters used for the calculation of the quantum phase fluctuations are 
defined as 



U {e,t, = (AiV)2 [(A^)2 + (AC) 



(1.37) 



ay] ={/\N)\/\SY (1.38) 



and 



Q (e, t, lap) = S {e, t, lap) /< C >2 (1.39) 

where S and C are sine and cosine operators respectively (explicit definitions are given in chapter 
3). Reduction of these quantum phase fluctuation parameters from their initial values is possible 
for some particular values of 9 and t. When S{9, t, |ap) gets reduced to below its coherent state 
value then at least one nonclassical effect (either magnetically squeezed electromagnetic field or 
antibunching of photons or both) is observed. 

Reduced fluctuation states can not be represented by P-representation. Therefore, they are in 
the paradigm of nonclassical states. In the present thesis we study the possibilities of observing 
these nonclassical states in the interaction of a coherent electromagnetic field with a nonlinear 
medium. From the above examples it is clear that in order to study the possibilities of observ- 
ing nonclassical states we have to calculate quantum fluctuations in different observables (e.g. 
N, S, C, X,X etc.). Again if we have to calculate quantum fluctuations in any of these observ- 
ables then we have to know the time evolution of creation and annihilation operators. Therefore, 
we have to obtain operator solution of the equation of motion of the anharmonic oscillator. This 
is why, we have to work either in Heisenberg picture or interaction picture. In the next chapter 
we obtain the operator solutions and use them to study different quantum optical and nonlinear 
optical phenomena in the subsequent chapters. 

Chapter 2 deals with the operator solution of the equations of motion ( |1.30D corresponding to 



our model Hamiltonian ( |1.29| ). The works reported in this chapter are arranged according to the 
time of their appearance to give a clear exposure of the time development of the subject as well 
as that of our work. For example, at first we report classical and quantum solutions of quartic, 
sextic and octic anharmonic oscillators by using Taylor series technique. While working with 
Taylor series, we observe that the approach is lengthy and time consuming. But many intrinsic 
symmetries of the problem is made transparent to us. For example, we observe that first order 
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correction to the frequency operator is always a function of the unperturbed Hamiltonian (-ffo)[|- 
Using these observations and imposing a physical condition that a correct solution should give 
correct frequency shift we have succeeded to generalize the results of Bender and Bettencourt 
[ pO| ] and have obtained the solution of the generalized AHO [|37|]. Later on we have general- 
ized the quartic oscillator solutions obtained by various other techniques and observed that an 
apparent discrepancy is present between the solutions obtained by different techniques. Then the 
question arises: Which solution should be used for our physical calculations? Therefore, at the 
end of chapter 2 we compare different solutions and conclude that all the correct solutions are 
equivalent and the apparent discrepancy is due to the use of different ordering of the operators. 
Showing this equivalence we become mathematically equipped to investigate the possibilities of 
observing different nonlinear optical phenomena in a nonlinear dielectric medium. 

In chapter 3 we study quantum phase fluctuations of coherent light. We begin by asking the 
question: How can one write down a quantum mechanical operator corresponding to the phase 
of a harmonic oscillator or equivalently a single mode of electromagnetic field? This question is 
the statement of the quantum phase problem. Since Dirac had introduced this problem in 1927 
[^2|l many people tried to provide a satisfactory Hermitian phase operator. As a result of these 
attempts there exist different definitions of the phase operator. In chapter 3 we shortly review 
a few of them and choose Pegg Bamett approach |44l ] to calculate the quantum phase fluc- 
tuation parameters S{9,t, \a\'^),Q{9,t, and U{9,t, by using the operator solutions of 
quartic anharmonic oscillator derived in chapter2. These fluctuation parameters were already 
studied by Gerry and Lynch for RW approximated Hamiltonian ( |1.33| ) of the same sys- 
tem. Here we observe that enhancement as well as reduction of S{9,t, Q{9,t, and 
U{9, t, with |a;p (number of photons present before the interaction) is possible for different 
values of the free evolution time t and phase 9 of the input coherent light. This observation is in 
sharp contrast with the earlier results |^]. Thus the importance of inclusion of photon number 
nonconserving terms in the model Hamiltonian is established. This interesting observation of 
chapter 3 has inspired us to study the other nonlinear optical effects with the same model. 

In chapter 4 we study the possibilities of generating electromagnetically squeezed light in a 
third order nonlinear medium. At first an analytic expression for the quantum fluctuation in phase 
quadrature X is obtained (up to the second power in A). We observe that an electrically squeezed 
electromagnetic field can be produced due to the interaction of an intense beam of coherent light 
with a third order nonlinear medium. A few special cases are discussed. 

In chapter 5 quantum statistical properties of the radiation field is discussed. To be more 
precise, bunching and antibunching of photons and the photon number distribution (PND) is 
studied. We observe that we can obtain nonclassical phenomenon of antibunching due to the 
interaction of an electromagnetic field with third order nonlinear medium of inversion symmetry 
for particular values of interaction time and phase of the input coherent state. 

In chapter 6 our studies on third order nonlinear medium is extended to the (m — l)-th order 
nonlinear medium in general. We study the nonadiabatic geometric phase or Aharonov Anandan 
phase for our physical system. Our works of chapter 4 and 5 are extended to the case of (m — 1)- 
th order nonlinear medium. Possibilities of generating nonclassical states in an (m — l)-th order 



An extensive proof of this observation is given by Speliotopoulos [MvL 
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nonlinear medium are also discussed. 

Finally in chapter 7 an outlook on the present work is given. In particular we talk about the 
limitations of our endeavor and future scope of investigations along the line of thought indicated 
by us. 
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Chapter 2 

Mathematical preliminaries 



Physical scientists have always tried to understand the Nature in terms of simple models. The 
simple harmonic oscillator (SHO) is perhaps the most useful one among them. A particle sub- 
ject to a restoring force proportional to its displacement gives rise to the model of a SHO. The 
Hamiltonian corresponding to a classical SHO of unit mass and unit frequency is given by 

H = ^- + '-. (2.1) 
2 2 ^ ' 

The Hamilton's equations for the SHO are 

^ " a/F"""" (2.2) 
X = ^ = p . 

Therefore, the equation of motion of the SHO is 

x + x = Q (2.3) 

where equation ( |2.2| ) is used. The solution of the equation ( [2. 3D is found to be 



xo(t) = x(0) cost + x(0) sin t. (2.4) 

The parameters x{0) and x{0) are the initial position and momentum of the oscillator. The 
subscript '0' denotes the zeroth order (i.e A = 0) solution. Thus the position of the oscillator at 
a later time t is completely known in terms of the initial position and momentum. Nevertheless, 



the momentum is also known from the Hamilton's equations (P3- Hence, the oscillator problem 
is completely solved. 

The equation of motion of a quantum SHO may simply be obtained by imposing x and x 
as noncommuting operators and the solution of a quantum SHO in Heisenberg picture is the 
operator equivalent of the equation ([2.4D. 

Now for the real physical problems, anharmonicity and/or damping are to be incorporated in 
the model Hamiltonian and hence in the equation of motion. However, in the previous chapter we 
have seen that only anharmonic term appears in the model Hamiltonian of our physical system. 
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The general form of the Hamiltonian of a classical anharmonic oscillator having unit mass and 
unit frequency is given by 



2 2m 



(2.5) 



where m (> 3) is an integer and A is the anharmonic constant. Depending upon the problem of 
physical interest, different types of anharmonic oscillators appear. Of course, the anharmonic 
constant A is different for different types of anharmonic oscillators. The equation of motion 



corresponding to the Hamiltonian ( [2.5[ ) is 



X + X + Xx 



m— 1 



0. 



(2.6) 



where the equation is used. We need an operator solution of ( |1.30| ) to study the quantum 
fluctuations of coherent light in nonlinear media. But equation of motion ( ]1.30D and its classical 
counterpart ( |2.6D are not exactly solvable in a closed analytical formfj. Keeping this in view, 
present chapter is devoted to the study of classical and quantum anharmonic oscillator problem. 
This problem is already investigated by several authors. To begin with, we recall the problem of 
a classical quartic oscillator (Duffing oscillator). 



2.1 Classical and quantum solutions of the quartic oscillator 



For m = 4, equation ( |2.6D has an exact solution in terms of the elliptic function [p3|]. This 
solution is available only in the phase plane. For positive anharmonic constant (A > 0), the 
solution is periodic with a fixed center. The period is also obtained in terms of the elliptic func- 
tions. However, this approach is not useful in many occasions. For example, it is not possible to 
predict the trajectory of the particle executing quartic anharmonic motion. For this reason, the 
solution in the phase plane is hardly useful except in the nonlinear dynamical studies. Keeping 
the above facts in view, a large number of approximate methods have been devised for the pur- 
pose of getting an analytical solution of the classical Duffing oscillator problem. These include 



perturbation technique [[4-5D, variation of parameters 0460 and Taylor series approaches [|47|]. The 
ordinary perturbation technique, a pedestrian approach, leads to the unwanted secular terms p5| ] . 
The removal of secular terms from the solution is a serious problem. There are some methods 
which are successfully used to sum up the secular terms for all orders. These include "tucking in 
technique" [Q, multiscale perturbation theory (MSPT) and renormalization technique [|45|]. 

The noncommuting nature of the operators poses serious difficulties for the purpose of getting 
approximate operator solutions to the quantum anharmonic oscillators. However, few methods 
are available for the purpose. We have already mentioned these methods pB| - PT| ] in 1.2.2. In the 
next two subsections of the present thesis we derive solutions of classical and quantum quartic 
anharmonic oscillators respectively. Actually, at first we develop a second order solution for 
the QQAHO using the Taylor series approach and after that, we construct a first order solution 



' Exact solution of ( [2.6[ ) can be obtained in terms of elliptic functions. However, that solution can not be used 
to predict the trajectory of the particle executing anharmonic motion. So the solution obtained in terms of elliptic 
functions are not exact in a strict sense. 
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for m — th oscillator in general using various other techniques. Solutions obtained by different 
techniques are compared in the subsequent subsections. The solutions are used to study the 
interaction of the strong electromagnetic field with a dielectric medium in the later part of the 
thesis. 



2.1.1 Classical and quantum solution of quartic oscillator 

Let us attempt a second order solution of the QQAHO by using the Taylor series approach de- 
vised by Mandal . At first we develop the solution of the QAHO in classical picture and after 
that a straight forward generalization is made for getting the corresponding solution in quantum 
picture. 

The differential equation of the classical Duffing oscillator of unit mass and unit frequency is 

X + x + Xx^ = (2.7) 
where the parameter A is the anharmonic constant of the system. The solution of the equation 



( P/TD may be written as the sum of different orders of anharmonicities. Thus the corresponding 
solution is 

x{t) = xoit) + xi{t) + X2it) + ... (2.8) 

where Xo{t), xi{t) and X2{t) are zeroth, first and second order solutions respectively. The first 
and second order solutions are proportional to A and A^ respectively. The zeroth and first order 
solutions are available from the works of Mandal [ p2| ] and are given by 

Xoit) = x{0) cost + i;(0) sint (2.9) 



and ^ 

Xi{t) = — '^^32°'' (cos t — COS 3t + 12t sin t) 



+ ^^^^§^(sin3t-7sint + 4tcost) 
22 (cos 3t — cost + 4t sin t) 
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(2.10) 



i^(0). 



32 (sin 3t + 9 sint — 12t cost) 

respectively. The purpose of the present subsection is to find the second order solution X2(t) of 
the classical quartic oscillator and to introduce the Taylor series method which is also used here 
to deal with the higher anharmonic oscillators. The solution of the classical Duffing oscillator 
(without forcing term) is given by the following Taylor series 

t^ 

x{t) = x{0) + tx{0) + —x{0) + . . . . (2.11) 

The assumed Taylor series is expanded about the origin. Furthermore, we assume that t is suf- 
ficiently small and is such that the above series ( |2.1 1| ) expansion is possible. Now we express 
higher order time derivatives of a;(t) at t = 

x-(0) = -x(0) -3Ax2(0)i;(0) 

x-(0) = a;(0)+4Ax3(0)-6Ax(0)i;2(0) + 3A2x5(0) 

x-(0) = x{0) + 24Aa;2(0)i;(0) - 6Xx^{0) + 27AV(0)i;(0) (2.12) 
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where we neglect the terms beyond and always substitute x by —x — \x^. 
The relations ( ]2.12D are substituted back in equation ( |2.1 ID and the coefficients of 



a;^(0), a;^(0)a;(0), a;3(0)i;2(o), a;2(0)i;3(0), x(0)x^(0) and x5(0) are collected. The correspond- 
ing coefficients are 

Ci = A2(3 - 51 + 846 - 15078 + 300705 - 6633081 + .. .) 

C2 = A2(27- 639 + 13230 - 284094 + ...) 

C3 = A2(126 - 3888 + 98820 - 2440224 + ...) 

= A2(378 - 14688 + 429084- 11416464+...) (2.13) 

= A2(756 - 33156 + 1023948 - 27952668 + .. .) 

and 

Ce = A2(756 - 33156 + 1023948 - 27952668 + .. .) 

where the contributions from Taylor expansion part are not taken into account. These contribu- 
tions are taken care of in the later part of this work. Constructing difference equations for the 
above series and solving them we can obtain the r — th term of the coefficient For example, 
the r — th term of Ci is given by 



t; = (-1)'+^ X ^ (25''+i + 216r x 9*" + 288r2 + 240r - 25) 



(2.14) 



Now we can write the r — th term of the coefficient of x (0) as 

^r-^r^ (2r + 2)!' ^^-^^^ 

where the factor comes from the Taylor expansion part. The net coefficient of a;^(0) is 

obtained as 

Ki = Y,'iLoU = ]^(cos5t — 36t sin3t — 24 cos3t n ^f^\ 

- 72t2cost + 96tsint + 23 cost). 

The parameter Ki exhibits fifth and third harmonic generations. The remaining coefficients are 
also obtained by using the similar procedure as adopted for the evaluation of Ki. The corre- 
sponding coefficients are 

K2 = Yi^(5sin5t+108tcos3t-132sin3t-72t2sint + 599sint-336tcost) 

= ^(-5cos5t + 90cos3t + 36tsin3t-72t2cost-85cost + 264tsint) 
K4 = ^(-5sin5t + 36tcos3t + 6sin3t-72t2sint + 427sint-456tcost) 
K5 = Y^(5 cos 5t + 108t sin 3t+ 108 cos 3t- 72)^2 cost -113 cos t + 240t sin t) 
and 

Kg = Y^(sin5t-36tcos3t + 48sin3t-72it2sinit + 271sint- 384tcost). 

(2.17) 

Hence the desired second order solution is given by 

X2{t) = KiX^iO) + K2x\0)x{0) + K3X^{0)x\0) 

+ Kix\0)x^{0) + K5x{0)x\0) + Kex^{0). ^ ^ 
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The total solution Q) is simply the sum of (jO^) and ( P^ . For i;(0) 
the total solution reduces to an extraordinarily simple form 



and x(0) = A, 



x{t) 



Acost - ^( 



COS t — COS 3t + 12t sin t) + (cos St 



36tsin3t - 24 cos3t - 72 cost + 96tsint + 23 cost). 



(2.19) 



It is clear that the presence of the secular terms proportional to t sin t, t sin 3t and cos t pose a 
serious difficulty as t grows. For weak coupling case {i.e A ^ 1), these terms may be summed 
up with the help of various summation techniques. However, for strong coupling the summation 
technique would fail and would give divergent solutions. We confine ourselves in the weak cou- 
pling regime and hence the summation for different orders of A is possible. For weak coupling, 
the secular term in the first order solution is known to produce the frequency shift of the oscilla- 
tor . In the present case similar frequency shift may be observed in following way. Equation 



(2.19) can be rearranged as 



x{t) 



A 



1 - 



2! 



cos3t- |AA2tsin3t) 



^8 
1024 



cost - A(^\AH) siYit-^ (cost - fAA^tsint) 



[cos 5t - 24 cos 3t + 23 cos t] + ^^^t sin t. 



(2.20) 

Now we use the tucking in technique to remove the secular terms. As we are dealing with the 
second order solution we can write 



sin(6At) 

cos(6At) 
{h'X) sin(6At) 
{h'X) cos(6At) 



6At 

_ {bXtf 
2! 

h'hXH 
h'X 



(2.21) 



where h and h' are constants and the terms beyond A^ are neglected. Thus the secular terms can 



be removed from the solution and finally the solution ( P.20D reduces to 



XA^ 



X^A'> 



x{t) = Acosu^t — — — (cosLU^t— cos3u;^t)+ (cos5u;^t— 24cos3a;^t+23cosu;^t). (2.22) 
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The shifted frequency is given by 



1024 



-XA^ - —X^A^ 
8 256 



(2.23) 



which coincides exactly with the frequency obtained by using Laplace transform [|49D. Thus, 
the existing second order solution of a QAHO p9| ] is reproduced as a special case (z.ex(O) = 
0, x(0) = A) of our solution (12. SP. For x(0) = 0,x(0) = A the frequency shift of the oscil- 
lator due to quartic anharmonicity is (^| — ^^)3A. The frequency shift due to the first term 
(SXA^/S) dominates over the frequency shift due to the second term (21A^A^/256) as long as A 
is small. Hence, the shift in the frequency increases with the increase of A till the contribution 
of the second term compensates the contribution of the first one. Further increase of A causes 
a frequency shift in the opposite direction. Thus the second order correction to the frequency 
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shift is quite consistent with the physical system. Now to have an idea of how good the solutions 
( |2.19D and ( [2.22D are we compare them with the exact numerical solution obtained by Mathe- 
matica (figure 1). From figure 1, we observe that the first order frequency renormalized solution 
( |2.22D coincides exactly with the exact numerical solution obtained by using Mathematica (for 
a = 2 and A = 0.05). However, the solution ( |2.19| ) diverges with the increase of time (figure 1). 
The divergent nature of the solution is a manifestation of the presence of secular term. 




2.1.2 Solution of quantum quartic oscillator 



It is straightforward to extend the classical solution ( |2.8D into a quantum solution for the anhar- 
monic oscillator by using the condition of quantization [ ]50| ] . The corresponding solution is given 
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by 



X{t) = X(0) cost + X(0) sint- ^^^(cost-cos3t+ 12tsint) 

+ 4[X2(0)X(0) +X(0)X(0)X(0) + X(0)X2(0)] X (sinSt- Tsint + 4tcost) 



^ [X(0)X2(0) +X(0)X(0)X(0) + X2(0)X(0)] X {cosSt- cost + 4t sin t) 
32 ^sinSt + 9sint — 12tcost) 
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+ ^-^^ (cos 5t - 36t sin 3t - 24 cos 3t - 72^ cos t + 96t sin t + 23 cos t) 
+ ^[X^(0)X(0) + X3(0)X(0)X(0) + X2(0)X(0)X2(0) 
+ X(0)X(0)X3(0) +X(0)X4(0)] 

(5 sin 5t + 108t cos 3t - 132 sin 3t - 72i(;2 sin t + 599 sin t - 336t cos t) 



X 



+ 5f^[X3(0)X2(0) + X2(0)X2(0)X(0) + X2(0)X(0)X(0)X(0) + X(0)X2(0)X2(0) 

+ X(0)X(0)X(0)X(0)X(0) + X(0)X(0)X2(0)X(0) + X2(0)X3(0) 

+ X(0)X(0)X(0)X2(0) + X(0)X2(0)X(0)X(0) + X(0)X3(0)X(0)] 

X (-5 cos 5t + 90 cos 3t + 36t sin 3t - 12t^ cos t - 85 cos t + 264t sin t) 

+ 5go[X3(0)X2(0) + X2(0)X2(0)X(0) + X2(0)X(0)X(0)X(0) + X(0)X2(0)X2(0) 

+ X(0)X(0)X(0)X(0)X(0) + X(0)X(0)X2(0)X(0) + X2(0)X3(0) 

+ X(0)X(0)X(0)X2(0) + X(0)X2(o)X(0)X(0) + X(0)X3(0)X(0)] 

X (-5 sin 5t + 36t cos 3t + 6 sin 3t - 12t^ sin t + 427 sin t - 456t cos t) 

+ 5f|y[X(0)X4(0) + X(0)X(0)X3(0) + X2(o)X(0)X2(0) + X3(o)X(0)X(0) + X4(0)X(0)] 
X (5cos5t + lost sin 3t + 108 cos 3t - 12t^ cost - 113 cost + 240tsint) 
+ (sin 5t - 36t cos 3t + 48 sin 3t - 72^^ sin t + 271 sin t - 384t cos t) 

(2.24) 

where X(0) and X(0) are the initial position and momentum operators. Of course, the direct 
use of the Taylor series for the equivalent operator differential equation ( |1.30| ) leads to the same 
solution ( |2.24| ). The quantization condition, [X(t), X(t)] = i, takes care the passage of classical 



28 



solution ( |2.8| ) to quantum solution ( |2.24D which may be written in the symmetrical form as 



X{t) = X(0) cost + X(0) sin t- 
+ §[X2(0)X(0)+X(0)X^ 



^^1^ (cos t — COS 3t + 12t sin t) 
X (sin 3t — 7 sin t + 4t cos t) 



+ 
+ 

X 

+ 

X 

+ 

X 

+ 

X 



§[X(0)X2(0) + X\0)X{0)] X (cosSt - cost + 4tsint) 

AX3(0) / 
32 V 
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sin 3t + 9 sin t — 12t cos t) 

(cos 5t - 36t sin 3t - 24 cos 3t - 12t^ cos t + 96t sin t + 23 cos t) 



^[X^(0)X(0)+X(0)X^(0)] 
(5 sin 5t + 108t cos 3t - 132 sin 3t - 72t^ sin t + 599 sin t - 336t cos t) 
|^[X3(0)X2(0) + X\0)X\0) + 3X(0)] 

(-5 cos 5t + 90 cos 3t + 36t sin 3t - 72t^ cos t - 85 cos t + 264t sin t) 
|^[X3(0)X2(0) + X2(0)X3(0) + 3X(0)] 

( -5 sin 5t + 36t cos 3t + 6 sin 3t - 72^^ sin t + 427 sin t - 456t cos t) 
^[X(0)X^(0)+X^(0)X(0)] 
(5 cos 5t + 108t sin 3t + 108 cos 3t 

A2X-'"'(0) 
1024 



fsin ht 



72t^ cos t - 1 13 cos t + 240t sin t) 
36t cos 3t + 48 sin 3t - 72^^ sin t + 271 sin t - 384t cos t) . 



(2.25) 

The equation ( [2 .251 ) is our desired solution for a quantum quartic anharmonic oscillator. The 
solution contains terms proportional to X'^ and X^ etcetera. In the classical limit these terms 
produce the cubic and fifth powers of amplitude respectively. Hence, as a limiting situation, 
the solution ( [2.25D gives rise to the solution corresponding to the classical QAHO. It is not at 
all surprising since the solution ( |2.25[ ) is obtained from its classical counterpart ( ]2.8| ) by the 
imposition of the quantization condition. 

The unpleasant secular terms can be removed (summed) from the solution ( |2.25D by using the 
tucking in technique. All the secular terms of equation ( |2.25| ) are tucked in to have the following 
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solution 

Xit) = J {x(0)cosvl>t + cosv[/tX(0)+X(0)sinvl/t + sinvl>tX(0) 

- ^^y^(cos^t-cos3^t) - ^(cos^t -cos3^t)X3(0) 

+ i[^^(0)^(0) +X(0)X2(0)] X (sin3^t-7sin^t) + §(sin3^t-7sin^t) 

X [X2(0)X(0) +X(0)X2(0)] - g[X(0)X2(0) +X2(0)X(0)] X (cos3*t-cos^t) 

- g(cos 3^t - cos m) X [X(0)X2(0) + X\0)X{0)] 

- ^^^^sin 3^t + 9 sin m) - 4 (sin 3m + 9 sin ^t)X3(0) 



+ ^^^24 (cos5^t-24cos3^t + 23cos^t) 

+ 1^24 (cos 5^t - 24 cos Sm + 23 cos m)X^{0) 

+ 2^[X4(0)X(0) +X(0)X4(0)] X (5sin5^t- 132sin3^t + 599sin^t) 

+ 2^48 (5 sin 5m - 132 sin Sm + 599 sin m) x [X4(0)X(0) + X(0)X4(0)] 

+ ^[X3(0)X2(0) +X2(0)X3(0) + 3X(0)] X (-5cos5^t + 90cos3^t-85cos^t) 

+ ^(-5cos5^t + 90cos3^t-85cos^t) x [X3(0)X2(0) + X2(0)X3(0) + 3X(0)] 

+ ^[^^(0)^^(0) +^^(0)^^(0)] X (-5sin5^t + 6sin3^t + 427sin^t) 

+ i&l-^ sin 5m + 6 sin 3m + 427 sin m) x [X^{0)X^{0) + X'^{0)X^{0)] 

+ ^[3X(0)] X (-5sin5^t + 6sin3*t + 403sin^t) 

+ ^(-5 sin 5m + 6 sin 3m + 403 sin m) x [3X(0)] 

+ 2^48 [^(O)-^^(O) +^^(0)^(0)] ^ (5cos5^t + 108cos3^t- 113cos^t) 
+ 2^(5cos5^t + 108cos3^t- 113cos^t) x [X(0)X4(0) + X^(0)X(0)] 
+ ^^|^(sin5*t + 48sin3^t + 271sin^t) 
+ Y^(sin 5m + 48 sin 3m + 271 sin m)X^{0) } 

(2.26) 

where n is the eigenvalue of the number operator (a^a). The parameter n may be regarded as 
the number of photons present in the radiation field corresponding to the electric field operator 
X(t). The prefactor of the above equation ( p.26D takes care of the proper normalization of X{t). 
The operator ^ is second order frequency operator and is given by 

^ = 1 + ^ - ^ {Q9H^o - 12X^(0) + I) (2.27) 

where Hq = — ^ H ^ is the unperturbed Hamiltonian. Now, the diagonal element of \E' in 

the number state basis is 

3A 1 

ijn =< n|^|n >= 1 + — (n + -) - —(51^^ + 5ln + 21) (2.28) 

4 2 o4 

where the following relations are used 

Ho\n > = (n + |) |n > 

X(0)|n> = ^ ^/n\n - 1 > +y/rrTl\n + 1 > forn^O' 

The consequences of our solution ( [2.26D may be compared critically with the existing results. 
For example, the energy and hence the frequency shift of a QQAHO is already known from the 



(2.29) 
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second order perturbation theory. The energy of a quartic oscillator inn — th state is given by 

3A 



[51] 



Er, 



1 3A 1 

fn + -) H (n^ + n + -) (34n^ + 51n^ + 59n + 21). 

^ 2^ 8 ^ 2' 128^ ' 



Therefore, the energy difference between two consecutive levels is 



AE = Er,- E, 



n-l 



1 + 



3An A^ 



64 



(51^2 + 21) 



(2.30) 



(2.31) 



In order to calculate the dipole moment matrix element < n— l\X{t)\n > and hence the shifted 
frequency of the quantum oscillator, we use equations ( [2.26[ ) and ( [2.28[ ). The corresponding 
calculated frequency shift coincides exactly with that of the frequency shift obtained by the 
perturbation method ( [2.3 1] ). The frequency shift for the vacuum field (i.en = 0) is evidently 
clear. It actually arises due to the vacuum field interaction. Here we can note that the Lamb shift 
is also explained in terms of the self interacting field (for a detailed discussion see reference [[52[]). 
Actually, our model Hamiltonian permits such an interaction through the quartic anharmonicity. 
Here we can note that a first order calculation is unable to predict the frequency shift for a vacuum 
field. 



2.1.3 Remarks on the solutions 

We obtain an analytical second order solution for a classical quartic anharmonic oscillator by 
using the Taylor series method. The solution is found to agree with the existing solutions obtained 
by various methods. The frequency shift of a QAHO increases linearly with the increase of 
anharmonic constant as long as the first order solution is considered [ |32| , p3[ |. These type of 
monotonic increase of the frequency shifts seem to be inconsistent. We obtain a correction to 
such monotonic increase of frequency shift. The solutions presented here are perturbative. We 
prefer perturbative solution instead of an exact numerical one. Because they can provide more 
physical insight into the problem. 

The classical solution is quantized in order to obtain the solution corresponding to a quantum 
quartic oscillator. Interestingly, the frequency shift of the quantum oscillator coincides with that 
of the frequency shift obtained by using a second order perturbation theory. The second order 
solution exhibits Lamb shift which is totally absent in a first order solution. Hence, we observe 
that our solutions reproduce the existing results for classical and quantum oscillators of quartic 
anharmonicity. 

From the mathematical point of view, the solutions presented in this section are of great aca- 
demic interests. In addition to these, the classical solution finds several applications in nonlinear 
mechanics. The solution corresponding to the quantum oscillator is useful in quantum optics and 
in the field theory. For example, the solutions are used to study the quantum statistical properties 
(e.g. squeezing, higher order squeezing, bunching and antibunching of photons) of radiation field 
interacting with a cubic nonlinear media in the subsequent chapters of the present thesis. More- 
over, the quantum phase of the output field is also studied with the knowledge of field operators 
X{t), X{t), a{t) anda"f(t). 
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2.2 Higher anharmonic oscillators : sextic and octic oscilla- 
tors 



In addition to the quartic oscillator, people have studied the higher anharmonic oscillators also 
[ p2| , p3| , |54t |. These studies are made to obtain the eigenvalues as function of anharmonicity. 
However, the operator solution of higher anharmonic oscillators were not reported since recent 
past. Recently we have reported first order operator solutions [ p9| ] of sextic and octic anharmonic 
oscillators by using the Taylor series method. We have also developed a first order solution 
for the generalized anharmonic oscillator p7|]. Fernandez obtained first order correction to the 
frequency operator for a few particular values of m []35|]. In the next two subsection of the 
present thesis, we present analytical solutions of classical and quantum oscillators of sextic and 
octic anharmonicities by using the Taylor series approach. In the later part of this thesis we 
discuss the generalized anharmonic oscillator and compare the results with that obtained in the 
Taylor series approach. The sextic and octic anharmonic oscillators have potential applications 
in the studies of nonlinear mechanics, molecular physics, quantum optics and in field theory. The 
solutions are used to obtain the frequency shifts of sextic and octic oscillators. The computed 
shifts are compared and found to have exact coincidence with the frequency shifts calculated by 
using a first order perturbation theory. 



2.2.1 Classical and quantum solutions of the sextic oscillator 



For m = 6, the Hamiltonian ( |23| ) and the equation of motion ( [2.6D correspond to the case of 



a sextic anharmonic oscillator. The solution for such an oscillator is obtained as the sum of 
different orders of anharmonicities. The corresponding solution is 

x{t) = xoit) +xi(t) + . . . (2.32) 

where xo(t) and xi{t) are zeroth and first order solutions respectively. The zeroth order solution 



( [2.4] ) is simply obtained for A = in equation (|2.6|). The purpose of the present section is to 
find the first order solution xi(t). By using the Taylor series approach introduced in connection 
with the quartic oscillator we can obtain first order solution of the classical sextic anharmonic 



oscillator as 



where 



5 

xi{t) = Y,LiX^^~\0)i' (2.33) 

i=0 

Lq = g|j(cos5t + 15cos3t — 16cost — 120tsint) 

Li = Jj(5sin5t + 45sin3t-280sint + 120tcost) 

L2 = ^(-5cos5t - 15cos3t + 20cost - 120tsint) 

L3 = ^(-5sin5t+15sin3t-140sint + 120tcost) (2.34) 

L4 = Jj(5cos5t-45cos3t + 40cost- 120tsint) 

and 

L5 = ,^(sin5t - 15sin3t - SOsint + 120tcost). 



The total solution ( |2.8| ) is simply the sum of ( |2.4|) and ( |2.33| ) 
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Now to have a check of the solution ( |2.32| ), we consider a special case i(0) = and a;(0) = a. 
Hence the total solution reduces to an extraordinarily simple form as 



x(t) = a cost H (cosSt + 15 cos St — 16 cost — 120tsint). 



(2.35) 



It is clear that the presence of a secular term proportional to t sin t makes life difficult as t grows. 
We confine ourselves to the weak coupling regime and assume C is a constant such that 



sin (CXt) 
cos (CXt) 



CXt 
1 



(2.36) 



Equation ( |2.36| ) is actually a first order form of ( |2.21[ ). Now, the equation ( |2.35D can be rearranged 
as 

5 15 1 

x{t) = acos[t(l + —Xa'^)] + Aa^(-— cost + — — cosSt + — - cos5t) (2.37) 
lb 24 12o oo4 



where the equation ( [2.36[ ) is used. Thus the secular term present in the first order solution ( |2.35D 
is summed up for all orders and the frequency shift of the oscillator is obtained. The corre- 
sponding frequency shift of the sextic anharmonic oscillator is ^Aa^. The shifted frequency 
of the oscillator may be viewed as the renormalization of frequency by the anharmonic interac- 
tion. Interestingly, the frequency shift obtained in the solution ( |2.37| ) exactly coincides with that 



obtained by Dutt and Lakshmanan [p4|]. 

Now, a quantum solution of the sextic anharmonic oscillator may be obtained by replacing 
x{t) and x(t) by noncommuting operators X{t) and X{t) respectively and imposing the com- 
mutation relation ( |1.23p . Corresponding operator solution of the sextic oscillator is 



+ 
+ 

X 

+ 
+ 
+ 

X 

+ 
+ 
+ 

X 

+ 
+ 

X 

+ 



X{t) = X(0) cos t + X(0) sin t+ 3^X^(0) X (cos 5t + 15 cos 3t- 16 cost- 120t sin t) 
^[X4(0)X(0)+X3(0)X(0)X(0) 
X\0)X{0)X\0) + X(0)X(0)X3(0) + X{0)X\0) 
(5 sin 5t + 45 sin 3t - 280 sin t + 120t cos t) 
^ [X^{0)X\0) + X\0)X\0)X{0) + X2(0)X(0)X(0)X(0) 
X(0)X2(0)X2(0) +X(0)X(0)X(0)X(0)X(0) +X(0)X(0)X2(0)X(0) 
X2(0)X3(0) + X(0)X(0)X(0)X2(0) + X(0)X2(0)X(0)X(0) + X(0)X3(0)X(0) 
(-5 cos 5t - 15 cos 3t + 20 cos t - 120t sin t) 

[X3(0)X2(0) +X2(0)X2(0)X(0) +X2(o)X(0)X(0)X(0) +X(0)X2(0)X2(o) 

X(0)X(0)X(0)X(0)X(0) +X(0)X(0)X2(0)X(0) +X2(0)X3(0) 
X(0)X(0)X(0)X2(0) +X(0)X2(0)X(0)X(0) +X(0)X3(0)X(0)" 
(-5 sin 5t + 15 sin 3t - 140 sin t + 120t cos t) 

[X(0)X^(0) -fX(0)X(0)X3(0) +X2(0)X(0)X2(0) 
X3(0)X(0)X(0) + X^(0)X(0)" 
(5 cos 5t - 45 cos 3t + 40 cos t - 1 20t sin t) 
^f2^(sin5t- 15sin3t-80sint + 120tcost) 

(2.38) 
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where X{0) and X(0) are the initial position and momentum operators. The solution ( |2.38[ ) may 
be written in symmetrical form as 



X{t) = X(0) cost + X(0) sin t 

+ ^X5(0)(cost + 15cos3t - 16 cost - 120tsint) 

+ 4[X^(0)X(0) + X(0)X4(0)] 

X (5 sin 5t + 45 sin 3t - 280 sin t + 120t cos t) 

+ ^[X3(0)X2(0) + X2(0)X3(0) + 3X(0)] 

X (-5 cos 5t - 15 cos 3t + 20 cos t - 120t sin t) 



2\ 
768 



768 



[X\0)X\0) + X\0)X\0) + 3X(0)] 
-5 sin 5t + 15 sin 3t - 140 sin t + 120t cos t) 
[X(0)X4(0) + X4(0)X(0)] 



X (5 cos 5t - 45 cos 3t + 40 cos t - 120t sin t) 

+ ^^f68°^ (sin 5t- 15sin3t-80sint + 120tcost). 



(2.39) 



Equation ( [2.39[ ) is our desired solution for a quantum sextic anharmonic oscillator. Now we 
rearrange our derived solution (|2.39[) with the help of (|2.36D as 



X{t) 



+ 

+ 
+ 
+ 
+ 
+ 
+ 



i {X(0) COS (t + f {h^, + i}) + cos (t + f {h^, + i}) X(0) 
X(0) sin (t + ^ {h^, + 1 }) + sin (t + ^ {h', + i}) X(0) 

cos 5t + 15 cos 3t — 16 cos t) 
^ [X^(0)X(0) + X(0)X^(0)] (5 sin 5t + 45 sin 3t - 280 sin t) 



2\X''{0) 
384 

384 

^ [X3(0)X2(0) + X2(0)X3(0) + 3X(0)] (-5 cos 5t - 15 cos 3t + 20 cos t) 
^[X3(0)X2(0) +X2(0)X3(0) + 3X(0)]( 
3^[X(0)X^(0) + X^(0)X(0)](5 cos 5t - 45 cos 3t + 40 cost) 



-5 sin 5t + 15 sin 3t - 140 sin t) 



384 



[sin 5t — 15 sin 3t — 80 sin t) | 



(2.40) 

is the unperturbed Hamiltonian. Now the matrix element < n — 
l|X(t)|n > and hence the frequency shift of the oscillator may be calculated with the help of 



where H< 







+ 



X^(0) 
2 
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(|2.29D. Now, the dipole matrix element in terms of the number eigenket \n > is 



<n-l\X{t)\n> = ^2!(|^|<^_i|x(0)|n>cos(t + ^{n2 + i}) 

+ cos{t + ^[n^ + 1}) <n-l\X{0)\n> 

+ <n-l|X(0)|n>sin(t + ^{n2 + i}) 

+ sin(t+^{n2 + i}) <n-l|X(0)|n> 



+ <n 

A 



^'^^4°"' ("^^^ ~^ '^^^ 3t — 16 cos t) 



+ 3|4[X^(0)X(0) +X(0)X^(0)](5sin5t + 45sin3t- 280sint) 
+ ^[X\0)X\0) + X\0)X^{0) + 3X(0)](-5 cos5t - 15 cosSt + 20 cost) 
+ ^[X3(0)X2(0) +X2(0)X3(0) + 3X(0)](-5sm5t + ISsinSt - 140 sin t) 



384 



[X(0)X^(0) + X^(0)X(0)] (5 cos 5t - 45 cos 3t + 40 cos t) 



+ ^^§^(sin5t- 15sin3t-80sint)] |n >}. 

(2.41) 

Hence, we obtain the frequency shift of the oscillator as ^{n^ + |)- Note that the frequency 
shift of the oscillator is observed even when there are no photons (i.e vacuum field) present in 
the radiation field. The vacuum field interacts with the medium and causes the frequency shift. 
This phenomena is a direct outcome of pure quantum electrodynamic effect and has no classical 
analogue. The frequency shift of the oscillator for a vacuum field is actually a second order effect 
and is discussed in the context of the second order operator solution of the quartic oscillator. The 
solution ( |2.4UD is consistent with the generalized solutions obtained in the later part of this thesis. 



At present, the frequency shift in ( [2.41D may be compared with the frequency shift calculated by 
using the first order perturbation technique [jSl]] . As a check, we calculate the energy of a sextic 
oscillator in the n — th eigenstate which is given by 

1 5A 

En = {n + -) + — (4n=^ + Qn^ + 8ri + 3) . (2.42) 

2 4o 

The energy difference between two consecutive levels is 

5A 1 

^E = En- En-1 = 1 + — + -). (2.43) 
Thus our calculated frequency shift and hence the solution agrees nicely. Finally we normalize 
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the equation ( [2 .401 ) with the help of the equation ( |2.41| ) as 



X(t) 



2cos('- 



W- {x(o)co. (t + f {hi + 1}) + CO. (t + f {hi + i}) X(0) 



+ 

+ 
+ 
+ 
+ 
+ 
+ 



X(0).m (t+'-f {hi + \])+sm{t+'-f {h^ + i}) X(0) 



2AX5(0) 



384 



cos5t + 15 cos'it — 16 cost) 



384 
2A 



[X^(0)X(0) + X(0)X4(0)](5sm5t + 45sm3t - 280smt) 



^[X3(0)X2(0) + X2(0)X3(0) + 3X(0)](-5cos5t - 15cos3t + 20cost) 



384 
2A 



^[X3(0)X2(0) + X2(0)X3(0) + 3X(0)](-5sm5t + 15sm3t - 140smt) 



384 

3^[X(0)X4(0) +X4(0)X(0)](5cos5t 

2AX5(0) , 
384 ' 



i5t — 15sm3t — 80smt)| . 



45cos3t + AQcost) 



(2.44) 

The equation ( [2.44] ) is the desired solution for quantum sextic anharmonic oscillator where the 
secular terms for all orders are summed up. 



2.2.2 Classical and quantum solutions of the octic oscillator 

Depending upon the nature of nonlinearity, the higher anharmonic oscillators come into the pic- 
ture. For example, the Hamiltonian and the equation of motion of classical octic anharmonic 
oscillator may be obtained if we put m = 8 in the equations ( |2.5| ) and ( |2.6D respectively. Using 
the similar procedure (as adopted for sextic oscillator), the solution (up to the linear power of A) 
for classical octic oscillator follows immediately as 



x{t) = 



+ 
+ 
+ 
+ 
+ 
+ 
+ 
+ 



x{Qi)cost + x{Q)sint 

^g^(cos7t + Ucosht + 126cos3t - 141cost - MQtsint) 
^^^^(7sm7t + 70sm5t + 378smt - 2373smt + Wdtcost) 
^^^^^(-21cos7t - 126cos5t - 126cos3t + 273cost - 2520tsint) 
^^^^^(-35sm7t - 70szn5t + 630sm3t - 3815smt + 2520tcost) 
^2-^^(35cos7t - 70cos5t - 630cos3t + mhcost - 2520tsint) 
^2-^^(21sin7t - 126sm5t + 126sm3t - 2415smt + 2520tcost) 
^^^f^(-7cos7t + 70cos5t - 378cos3t + 315cost - SAOtsint) 
^^(-sm7t + Usin5t - 126sin3t - 525sint + 8AQtcost) 



(2.45) 



To check the validity of our solution, we consider a special case x(0) = b and x(0) = 0. After a 
little algebra we have 

35A 

x{t) = bcos{l + j^b^)t + ^^(-141cost + 126cos3t + Ucos5t + cos7t). (2.46) 

Thus the frequency shift of the classical octic oscillator is proportional to the sixth power of 
amplitude as long as the first order solution is concerned. The solution (|2.46D as well as the 
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frequency shift have exact coincidence with the solution obtained by the procedure introduced by 
Bradbury and Brintzenhoff [|55|]. The solution of the classical octic oscillator ( |2.45D is now used 



to obtain the corresponding solution for its quantum counter part. It is found that the solution 



contains the secular terms similar to those appeared in ( [2.24D . The problem of secular terms are 
taken care by the same procedure as it is done in the case of a sextic anharmonic oscillator and 
the final solution appears as 

^(t) = o f35Al 2^^^.1 {x{0)cosQt + cosQtX{0) + X{0)sinQt + sinQtX{0) 

ZCOSl „ , ( DTi +0 It 1^ 



^^^(cos7t + Ucos5t + 126cos3t - Ulcost) 



^[X{0)X%0) + X^{0)X{0)]{7sin7t + 70sin5t + 378sint - 2373sint) 
^[X'{0)X\0) + X\0)X^{0) + 10X3(0)] 
(-21cos7t - 12Qcos5t - 12Qcos3t + 273cost) 
^[X\0)X%0) + X\0)X\0) + 9(X2(0)X(0) + X{0)X\0))] 
{—3bsin7t — 70sin5t + 630sm3t — 3815smt) 
^[X%0)X\0) + X\0)X%0) + 9{X\0)X{0) + X(0)X2(0))] 
{35cos7t — 70cos5t — 630cos3t + 665cost) 
3^[X5(0)X2(0) + X\0)X^{0) + 10X3(0)] 
(21sm7t - 12Qsin5t + 126sin3t - 2415sint) 

^[X%0)X{0) + X{0)X%0)]{-7cos7t + 70cos5t - 378cos3t + 315cost) 
+ ^'3072°'' {—sin7t + 14sm5t — 126sm3t — 525smt)| 



+ 

X 

+ 

X 

+ 

X 

+ 

X 

+ 



(2.47) 

where the first order frequency operator is ^7 = 1 + ^{AHl + hHo). Now, the matrix element 
< n — l\X{t)\n > is calculated to obtain the shifted frequency as cu' = l + ^{'n? + 2n). Interest- 
ingly, the frequency shift ^(n^ + 2n) depends on n^. It is clear that the vacuum field has no role 
to play in the frequency shift of the quantum octic oscillator. To have a possible verification of 
the solution ( |2.47| ), we calculate the frequency shift by using Rayleigh- Schrodinger perturbation 
theory. A first order calculation under this theory shows that the energy of the octic oscillator in 
the n — t/i state is given by 

En = {n + -) + — (- + 4n + Sn^ + 2n3 + n^). (2.48) 
2 d4 2 

The energy difference between two consecutive levels is /\E = En — -En-i = 1 + ^{n^ + 2n). 
Thus, the frequency shift coincides exactly with the frequency shift obtained by us. 



2.2.3 Remarks on the solutions: 



We obtain analytical solutions (first order) for sextic ( [2.33[ ) and octic ( [2.45| ) oscillators in classical 
pictures by using the Taylor series method. The solutions exhibit the presence of the usual secular 
terms. The secular terms for all orders are tucked in to obtain the frequency shifts of those 
oscillators. It is found that the calculated shifts agree satisfactorily with the shifts obtained by 
other methods. 
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The position and momentum of classical oscillators are replaced by their corresponding op- 
erators and the commutation relation is imposed. Hence, the solutions for quantum sextic and 
octic oscillators in the realm of the Heisenberg picture is obtained. Interestingly, the vacuum 
field has no role to play in the first order frequency shift of the quantum oscillators of quartic and 
octic anharmonicities. However, the vacuum field causes the frequency shift of a quantum sextic 
oscillator. 

Till now we have discussed the Taylor series method for obtaining first order operator so- 
lution of the quartic, sextic and octic anharmonic oscillators. A close look into the difference 
equations and their solutions developed in the earlier sections shows some symmetries among 
them. Those symmetries encouraged us to search methods of obtaining solutions of m — t/i 
anharmonic oscillator in general. In the next section of this chapter we develop a zeroth order 
MSPT solution for the m — th anharmonic oscillator and generalize the first order operator solu- 
tions of quartic oscillator obtained by various other techniques to obtain the first order operator 
solution of generalized anharmonic oscillator in the subsequent sections. One more thing we 
would like to note here, the Taylor series and MSPT methods for obtaining solutions of ( |1.30D 
for particular values of m are difficult and tedious. So in the next few sections of this chapter we 
shall try to develop simpler techniques for solving (|1.30[) in general. 



2.3 Generalized quantum anharmonic oscillator using an op- 
erator ordering approach 

In the present section we study the generalized quantum anharmonic oscillator problem in the 
Heisenberg representation by making use of two simple normal ordering theorems for the expan- 
sion of (a + a))"'. We observe that the merit of the present work is its simplicity. For example, 
the present method provides, on the one hand, a straight forward mathematical frame work to 
construct expressions for the energy eigenvalues as well as frequency shifts and, on the other 
hand, generalizes some of the recent results on the topic by Bender and Bettencourt [BD|-BT]]. In 



the subsequent sections of these chapter we have shown that the above mentioned theorems may 
be used to generalize the results obtained by other methods. 

We have already mentioned that the equation of motion corresponding to the Hamiltonian 



( |1 .29D of a generalized anharmonic oscillator having unit mass and unit frequency is given by 
( |1 .30D which can not be solved exactly for m > 2. However, a large number of approximation 
methods are available for solving ( |1.30| ) for particular values of m. We have discussed the Taylor 
series method in detail and have mentioned some other methods. MSPT is one of those meth- 



ods. Bender and Bettencourt pO|j3Tl ] generalized the existing theory of MSPT into an operator 



approach to get the zeroth order solution involving a quantum operator analogue of the classical 
first order frequency shift. Following the interpretation of Aks [ |28|j29| ] they interpreted it as an 



operator mass renormalization that expresses the first order shift of all energy levels. They ob- 
tained MSPT operator solution of a quantum quartic oscillator. In the present work we suggest 
a novel method in which MSPT results for anharmonic oscillators in general can be found. To 
that end we prove in the next subsection two theorems to construct a normal ordered expansion 



38 



of (a + a^)™ . In the subsequent subsections we obtain a generalized expression for the energy 
eigenvalues and a generalized solution for the equation of motion ( |1.30D . We also specialize 



our result to reproduce some of the existing results as useful checks on the generalized solution 
obtained by us. 

2.3.1 Operator ordering theorems 

On a very general ground one knows that in quantum mechanics and quantum field theory proper 
ordering of the operators plays a crucial role. Let /(a, a^) be an arbitrary operator function of 
the usual bosonic annihilation and creation operators a and a'^ which satisfy the commutation 
relation 

[a,a^] = l. (2.49) 

One can write /(a, a^) in such a way that all powers of always appear to the left of all powers 
of a. Then /(a, a^) is said to be normal ordered. In this work we want to write (a + a^)™ in 
the normal ordered form for integral values of m. Traditionally, for a given value of m this is 
achieved by using a very lengthy procedure which involves repeated application of ( |2.49| ). One 



of our objectives in this work is to construct a normal ordered expansion of (a + a^)™ without 
taking recourse to such repeated applications. 

We denote normal ordered form of f hy f^. On the other hand : / : denote an operator 
obtained from / by arranging all powers of to the left of all powers of a without making use of 
the commutation relation ([2.49D. Now if / = aa^ then = a^a + 1 and : / := a^a. Therefore, 



we can write 

: (a + at)™ :=: (a^ + ■= a™ Cia^a"^"' + ... +" aa^'a""' + •• + a^"^ . (2.50) 

Thus in this notation : (a + a'l")'" : is simply a binomial expansion in which powers of the are 
always kept to the left of the powers of the a. To write (a'^ + a)^ we can proceed by using the 
following theorems. 

2.3.1.1 Theorem 1. 

: {a^ + a)"" : {a^ + a) =: (a^ + a)"'+^ : +m : {a^ + a)™-^ : . (2.51) 
Proof: From ( |2.50| ) : (a^ + a)™ : (a^ + a) can be written in the form 

: (at + a)™ : (a^ + a) 

= at"""' + ('"Ci Co)at'"a + (™C2 Ci)a^'^~'a^ + .... + ('"C, Cr-i)a^"^~''^' a' + 
+ (™Ciat'""' +2"C2at'""'a + ... +r'™aat""''a'"^ + •••) 



m— 1 . 



(at + a)'"+i : +m : (a^ + a) 



In the above we have made use of the identities 



(2.52) 



(a''"at)jv = a^d'' + ra^' ^, 

r'^Cr = n'^-^C-i, (2.53) 

and "^+1^+1 = {""Cr Cr+l) . 
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Note that first identity in ( [2.53D can be proved with the help of the general operator ordering 
theorems [ p6| ] while the other two are trivial. 



2.3.1.2 Theorem 2 : 

For any integral values of m 



r=0,2,4.. 



t I ^\m-r 



with 



r! 
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(2.54) 



(2.55) 



(2.56) 



+ 1260 : (at + a)^ : +945 : (at + a) : . 
From (p.56D we venture to identify the general form of the above expansion for a given value of 



m as 



)N 



r=0,2,4.. 



'Cr : (at + ay 



(2.57) 



It is easy to check that ( |2.57[ ) gives all the expansions of ( [2.56D . So ( |2.57p is true for 



m 



1, 2, .., 9. To ensure the general validity of ( |2.57| ) we can use the method of induction. First we 



- F M Fernandez [ p^ has pointed out that this theorem can be proved by expanding both sides of the Campbell- 
Baker-Hausdroff formula [pq] 



about C = 0. 
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assume that ( |2.57| ) is true for a particular value of m. Now, 



(at + a)^+i 



(flt + a)^ (at + a) 



(at + a)"" : +"^€2 : 



a)™-2 : +... 



+ tr-2 ™C; 
= Z^r [tr-2 



2- [a 

r-2 ■■ lat + a)™-^+2 . m^^ . 

"'Cr-2{m-r + 2)+tr"'Cr] : 



(at + a)™-^ : +. 



at + a) 



at + a 



\m— r+l 



: (at + a) 



.m+l— r 



(2.58) 

Equation ( |2.58D exhibits that if ( |2.57[ ) holds for any arbitrary m it must hold for m + 1. This 
establishes that ( |2.57D gives the normal order expansion of (at + a)™ for any arbitrary integer m. 



We use this normal ordered expansion to study the anharmonic oscillator problem. 



2.3.2 Energy eigenvalues 

The first order energy eigenvalue Ei is < n\H\n >, where |ri > is the number state. In the 
literature there are two lengthy procedures [ ]5T| ] to obtain Ei. The first one is the usual normal 
ordering method. This method involves iterative use of ( [2.49D and the number of iterations in- 
creases very fast as m increases. For a given value of m we need [2™ — (m + 1)] iterations. Just 
to get a feeling of the number of iterations required for large values of m we can look at table 1 
below. 



m 


no. of iterations required 


2 


1 


6 


57 


10 


1013 


20 


1048555 


100 


1267650600228229401496703205275 



Table 2.1 

The table shows that for large m the construction of the expansion (at + a)^ becomes 
formidable. In the second procedure one proceeds by making repeated application of x operator 
on the number state. As in the normal ordering method this procedure is also equally lengthy. 
In the following we implement theorem 2 to derive an uncomplicated method to construct an 



expansion for Ei for the Hamiltonian ( |1.29D . 

Using the result in ( |2.54| ) we can write the expression for Ei for any integer value of m in the 
form 



El 



(n + i) + ^<n|E: 



r=0,2,4. 



U'^Cr : (at + a 



[n + ^} + 



r=0,2,4. 

J. m/^ m—r/^ 



m — r Qj 
2 



\n > 

a 2 \n > 



(2.59) 



m — r 
2 



m — r 
2 



' m—r ^ I 
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The expression ( |2.59| ) involves summations which are easy to evaluate and thereby avoids the 
difficulties associated with earlier iterative procedures. Although somewhat forced, it may be 
tempting to compare the simplicity sought in our approach with the use of logarithm in a nu- 
merical calculation or use of integral transforms in solving a partial differential equation. It is of 
interest to note that expression similar to ( |2.59[ ) can also be constructed for the first order energy 
eigenvalue of a Hamiltonian in which the anharmonic term is a polynomial in X. 



2.3.3 MSPT solution of the generalized quantum anharmonic oscillator 

We want to obtain MSPT operator solution of equation ( |1.30D for arbitrary integral values of m. 
The essential idea behind our approach is that the quantum operator analogue of the classical first 
order frequency shift is an operator function ( ^{Hq) ) of the unperturbed Hamiltonian Hq and 
secondly a correct solution should reproduce the first order energy spectrum. Now the general 
form of the zeroth order solution of ( |1.3C1| ) is given by 



Mt) = ^[X(0)cos(t + A(^(i/o)t) + cos(t + An(i/o)t)^(0) 
+ X{0)sin{t + Xn{Ho)t) + sin{t + Xn{Ho)t)X{0) 



(2.60) 



where G{n) is a normalization factor. So our task is to find out ^1{Hq) and G{n) in general. From 
equation ( |2.59D we obtain the energy difference for two consecutive energy levels as 

^n,n—l (-^l)m,n (-^l)m,n— 1 

= l + \uj{m,n) (2.61) 

1 I A ■^(■m— 2) , m/^ m—ri^ n—lr^ I m—r \ \ 

2^m 5 V 2 ' 

Since the correct quantum operator solution has to give the first order energy spectrum, we 
should have 

< n - l|Xo(t)|r;. > = < n - l|Xo(0)|n > cos[t + Acu(m, n)t] (2 62) 

+ < n — l\XQ{Q)\n > sin[t + \ijj{m^n)t\ . 

Equation ( |2.60D and ( |2.62| ) impose restrictions on our unknown functions Vl{Hq) and G{n). The 
condition imposed on Vl{Hq) is 

< n\n{HQ)\n > + <n - l\Q{HQ)\n - 1 >= 2uo{m,n) (2.63) 

or, 

+ -) + n(n - -) = 2cj(m,n) . (2.64) 

For a particular m the right hand side is a known polynomial in n and our job is simply to find 
out VL{n + i). We obtain this as 



E(-l)"-^a;(m,fc) 

.k=0 



+ (2.65) 

22m 
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Substituting the functional form of n{n + |) or Vl{Ho) in ( [2 .601 ) if we impose condition ( |2.62D 
we get 



G{n) 



2cos 



(2.66) 



The results in ( [2.65D and ( |2.66D when substituted in ( [2.60D solve the generalized quantum anhar- 
monic oscillator problem. 



2.3.4 Specific results and tlieir comparison witli tlie existing spectra: 

Although we have solved the quantum anharmonic oscillator in general, just now it is not possible 
to compare our solution directly with any other result since the present study happens to be 
the first operator solution of the generalized anharmonic oscillator. Later on we develop other 
techniques to solve the anharmonic oscillator in general and compare those result with the present 
one. Here we calculate some specific results from our general expressions and compare them 
with the existing results. 
For m = 4 we have, 



3n 3 
T + 8 



3 1. 
4(^+2' 



Therefore, 



In terms of (2. 



and G{n) 

the total solution for the quantum quartic anharmonic oscillator is 



4H0 



X{t)\ 



m=4 



2cos{^) 



X{0)cos[t + '^Ho] + cos[t + ^Ho]X{0) 



+ X{0)stn[t + ^Ho] + szn[t + ^Ho]X{0) 



(2.67) 



(2.68) 



(2.69) 



This coincides exactly with the solution given by Bender and Bettencourt [ pO| , |3Tl ] and also coin- 
cides with the Taylor series solution (|2.26D . This also gives the correct classical frequency shift 
[ p7| ] in the limit x(0) = a, x(0) = 0. Similarly we have 



X(t)\m=6 
x(t) \m-- 

and 

x{t)\m=W 



2cos(^n) 



XiO)cos[t + f + + cos[t + 'i{Hl + \)t]X{Q) 



+ X{Q)sin[t + 'i{Hl + \)t] + sin[t + f (i^^ + \)t]X{Q) 



(2.70) 



2cos 



^(6n2+3)t] 

X(0)cos[t + ^{AHl + 5iJo)t] + cos[t + ^-^{AHl + 5ifo)t]^(0) 
+ X{'d)sin[t + ^{AHl + 5/7o)t] + sin[t + ^{AHl + 5/7o)t]^(0) 



X 



+ cos[t + '-§{H^ + IHI + ^)t]X(0) + X(0).m[t + '^{H^, + iRl + ^)t] 
+ s^n[t+^^^^{H| + lHl + ^)t]X{Q) 



(2.71) 



43 



The solutions for sextic and octic oscillator exactly coincides with the solution obtained by us 
using Taylor series approach [|3^] and all solutions give the correct classical frequency shifts in 
the appropriate limit Q . 

2.3.5 Remarks on the solutions 

We conclude by noting that depending on the nature of nonlinearity in a physical problem the 
treatment of higher anharmonic oscillators assumes significance. But studies in such oscillators 
(for m > 4) are not undertaken in the Heisenberg approach presumably because the existing 
methods tends to introduce inordinate mathematical complications in a detailed study. In the 
present work we contemplate to circumvent them by proving a theorem for the expansion of 
(a^ + a)^. Thus the results of the present section are expected to serve a useful purpose for 
physicists working in nonlinear mechanics, molecular physics, quantum optics and quantum 
field theory. In fact, we use these solutions to study nonlinear and quantum optical effects in 
detail in the later part of the present thesis. 



2.4 Generalized quantum anharmonic oscillator using the 
time evolution operator approach 

Generalized quantum anharmonic oscillator can be studied through some alternative approaches. 
We have started our discussion on the possible alternative approaches from the MSPT approach. 
In the present section we study the time evolution operator approach. In this approach we find 
out the perturbed time evolution operator in the interaction picture Ui{t) as a function of a, a) 
and t and then we use that to find out the time evolution of the creation and annihilation operators 
which in turn solve the equation of motion ( 11.30| ) of the generalized anharmonic oscillator. The 
secular terms are removed from the solution by using the tucking in technique, discussed in the 
context of Taylor series solution. 

2.4.1 Time evolution operator in the interaction picture: 

In the interaction picture potential Vi and the time evolution operator f//(t) are respectively 

Vi{t) = exp{iHot)V exp{-iHot) (2.72) 

and 

Ui{t) = 1-1 [' dUViiU) + {-if f dhVj{ii) r dt2Vi{t2) + ... . (2.73) 

JO JO JO 

Where the suffix / stands for the interaction picture. Now to obtain the compact expressions for 



Vi{t) and Uj{t) (up to first order) for the model Hamiltonian (]6.16D of our interest we have to use 
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following operator ordering theorems 

Theorem 3: When / (a, can be expanded as a power series in a and then 



ex\i{qa)a)f {a^a)^ exp(— ^a^a) = / (^aexp(— <j), exp(<j)j (2.74) 

where ^ is a c-number [ ]5^ ] . 

With the help of theorem 2 (see section 2.3.1.2) and theorem 3 we can write 



Viit) = expUa'^at) — ^^(a'^ + a)"^ expi—ia'^at) 

m(2) "3" 



Er t2r "'Cs, Ei""^'^ (™-2r)^^^t^^™-2r-p g^p {tt{2p -771 + 2r)) 



(2.75) 



Now the first order time evolution operator in the interaction picture is 

Ui{t) = l-iJ^dt,Vj{h) 

= 1 - t—^ Jo (e? t2r ™C2. Ei^-^"^ 2-)Cpat''a— (iti(2j9 - m + 2r)) ) dh 

_ 1 A / Y^T J- mri x^{m-2r) lm-2r) ri ^ f ^ m-2r-p [cxp{it{2p-m+2r))- 1] \ 



m(2) 2 Y 2 



(2.76) 

The time evolution of the annihilation operator a in the Heisenberg picture an and in the inter- 
action picture a/ are related by 

o-uif) = Uj{t) exp{ia^at)a{0) exp{—ia^at)Ui{t) = exp{—it)ai(t). (2.77) 

Now the time evolution of the annihilation operator in the interaction picture is 

aj{t) = Uj{t)a{0)Ui{t) 



a- ' 



m(2)^ 



Xt 



Y^Tj- mr' Y^(™-2r) (m-2r)r^ ^fP-'^ ^m-2r-p p[cxp(it(2p-m+2r))-l] \ 



„ I 2 m-2r ^ mr^ {m~2r) „ n\ n — 9- 

~ ^ — TTST l-^r — o — '-2r *-'2r m-2r a' a 2 

m(2)^ V 2 

(2.78) 

So we have obtained the time evolution of the annihilation operator. 



2.4.2 First order frequency operator: 



In the solution of ( |1.30| ) frequency becomes a operator function. In the first order frequency 



operator is a function of the unperturbed Hamiltonian H^. The frequency operator is known for 
various values of m. Here we give a compact expression for the first order frequency operator 
Vl{Hq) and write the expression for some particular m to compare our result with the existing 
results. To do so we have to use the following identity 

a)"" = a'^ala'^a — l){a''a — 2)...{a^a — n + 1) 

r(go+j) (2.79) 
r(//o-n+i) • 
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The secular terms present in ( |2.78D can be tucked in into the expression of a by using the tucking 
in technique introduced in ( [2.21D . Now using ( |2.78| and |2.36D we can write 



a{t) = exp{—iXQt)a 

A I v^T + mr^ ^(m-2r) (m-2r)^ „t^-l exp{int{2p-m+2r-l)) oA\ 
~ ~^^^^\^'' '-^arZ.p.^m^ ' 'L^pCt' (2.80) 

^ ^m■-2r-p p[exp(^f(2p-m+2r))-l] ' 
2p— m+2r 



where 



(2.81) 



m(2)^ 2 y 

1 ( Y^T m-2rj jn/^ (m-2r)/^ r(ij'o+|) ^ 

m{2)^ V 2 J- (-Ho 2 ^2'/ 

Thus we have a compact expression for the first order correction to the frequency operator ( |2.81D 
and an expression for the time evolution of annihilation operator which is free from secular term 
( |2.80| ). Expression ( |2.8Q ) essentially solves the equation of motion ( ]1.3Q ) of the generalized 



anharmonic oscillator since 

X{t) = ^{a\t)+a{t)). (2.82) 

Here we must note that the frequency operator of the annihilation operator is not exactly the 
same with the operator usually used to write the first order frequency correction of the position 
operator. Relationship and equivalence criterion for the frequency operators obtained by using 
different techniques are discussed at the end of this chapter. 



2.5 Generalized quantum anharmonic oscillator using renor- 
malization group technique 

In this method a time parameter, additional to the initial value point is introduced, in such a way 
that the perturbation expansion is valid in the vicinity of the introduced time parameter. The 
coupling constants, constants of motion and/or initial conditions are suitably changed by the in- 
troduced time parameter. But the solution should not depend on the introduced time parameter so 
the first derivative of the perturbed solution with respect to introduced time parameter should be 
zero and this condition is known as RG condition. Solution of the differential equation imposed 
by RG condition can give provide us an operator solution of the anharmonic oscillator free from 
secular terms. 

From the perturbation expansion ( |2.78| ) the secular part of the annihilation operator up to the 



first order is 

,\{t — t) / ^ m — 2r 



m — 2r „ 
1 m — 2r 



-l- 



m(2)' 



^ ':^^t2rc2r ('""^'■^C^at^" a"V^ . (2.83) 



Now imposing the RG condition ^ = 0, we obtain 

( m 

^ m - 2r „ (m-2r)^ — \^^\ n 9.A\ 

> — t2r C'2r ' C m-2r a' a 2 . (2.84) 

^2 2 ' 
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Since d^a and [a,a^] are constants under the flow of r, we are aUowed to solve equation ( |2.84D 



as 



a[T) = exp —t- 



Xt 



^ m — 2r 



J. m/^ (m~2r) ^t" 
Zor '-^2r (_/ m-2r g 



m(2)f W 2 
Therefore the first order correction to frequency operator is 



a 2 



-1 a(0). (2.85) 



m{2)^ 



1 



m(2)"2" 



E2 m—2r+ mry (m~2r)ry 
r — — l2r *-'2r m-2r 



r(/fo+i) 



r(//o- 



t-2r I 3\ 
2 ^2-' 



(2.86) 



This is in exact accordance with the first order frequency operator derived by the time evolution 



operator method ( |2.81D 



2.6 Generalized quantum anharmonic oscillator using near- 
identity transformation technique 

The near-identity transformation relates the full solution for the annihilation operator a to the 
zeroth order term b as 

ait) = hit) + ATi(6t(t), hit)) + X^T^ih^t), hit)) + (2.87) 

and the equation for the time dependence of the zeroth order term in the normal form is 

^ = f/o(6t(t), hit)) + Af/i(6t(t), hit)) + X'U^ih^t), hit)) + .... (2.88) 



Now using the Heisenberg equation of motion for the Hamiltonian (|1.29| ) we obtain 



da 
dt 



ia + i4rria + a^)"'~^ 

m — 1 

2 



(2.89) 



Inserting equation ( |2.87| and |2.88D into ( |2.89D we obtain a relation between f/„ and T„ in each 
order n. These relations can not determine f/„ and T„ uniquely. To reduce this nonuniqueness 
problem we follow the work of Kahn and Zarmi and choose T„ not to depend on time 
explicitly. In the first order only the resonant terms of the form ih'^h)''h acts as the source terms 
for the explicitly time dependent part of the annihilation operator a. Thus from the expansion 



(2.89) we have 



1 



(f-i) 



2~ n 2 

r=0 



— 1 m—2r 



it)h^it). 



(2.90) 
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Substituting ( [2 .901 ) in ( [2.88[ ) we obtain the dynamical equation of b{t) as 



Therefore we have 

/ / A ^^^-^^ , , ^21^-1 m-2r A\ 

6 = exp I ( 1 + ^ X] ^2r '""'C2r ("^"^^-i) 6o ' j I (2.91) 

Now since the perturbation starts at t = so 6(0) = a(0) + C(A) and the zeroth order term h of 
the annihilation operator is 

h{t) = exp (^-it ^1 + ^ ^ E ^ ('"-^^-^)C^^2. ,at"^"Q^-i j j a(0). 

Therefore the first order correction to the frequency operator is 



O 1 2 -f m— 1/^ (m— 2r— 1)/^ „t" 



2^ 



a 2 



^Z^r=0 '^2r (--2r' ^ L- ^^^^^ m-ir J^j . (2.92) 

1 ( Y^T m—2r+ mr^ (m—2r)r^ „ r(_fi"o + 5) 



m(2)^ 



r -^t2r ^2t ^ -^^(u^Z^^r^) ) 



This is exactly equivalent with the frequency operator obtained by using the evolution operator 
technique ( |2.81| ) and the renormalization group technique (]2? 



2.7 Generalized quantum anharmonic oscillator using eigen- 
value approach 

This straightforward operator method is based on a generalization of the harmonic oscillator al- 
gebra for the creation and annihilation boson operators. This method was first used by Fernandez 
for the purpose of obtaining the first order correction to the frequency operator for particular val- 
ues of m. Later on we have generalized his idea to obtain a generalized expression for the first 
order frequency operator . The generalized treatment is given below. 

In this approach the Hamiltonian operator R is given and we look for an operator h such that 



[35-41] 



7 (2 93) 

Under such conditions it is possible to obtain eigenf unctions common to R and fi: 

= ^n^n. (2.94) 
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It follows from the equations ( |2.93| and |2.94| ) that 

{E^-E^ + nr,){^nWm)=^. 



(2.95) 



Since the operator h is nonzero, there must be a pair of states \E'm and '^n such that (\E'm|6|\E'„) 7^ 
0; therefore, 



Q = E — E 
The time evolution of the operator h is given by 

bit) = e'^^be-'^^ = ( 

therefore, 

(^nl&Wl^m) = e-**^"(vI/„|6|vI/„) = e^*(^"-^'-)(vI/„|6|vI/„). 
If we can write an operator O in terms of b and 6^ as O = 0{b, b'^) then we have 



(2.96) 

(2.97) 
(2.98) 

(2.99) 



If we rewrite equation ( |2.97[ ) as e**^^^^-'6 = 6e**^ and expand both sides in Taylor series 
about t = we obtain (H + fi)^6 = bH''. Consequently, for any operator function F(H) that 
can be expressed as an -power series we have [ ]4T| ] 

bF{H) = F{H + n)b, (2.100) 

which is consistent with equation (|2.96|). In order to generalize this result we can define 



b^'FiH) = F{k,H)b\ F{0,H) = F{H). 
A straightforward calculation shows that 

F{k, H) = F{k-1,H + n{H)). 
It is not difficult to prove that b'^b is a constant of the motion 

[H, b'^b] = b^VLb - b^VLb = 

but that 66^ is not 

[H, bb^ = -Qbb^ + bb^Q = [Q{H + Q{H) - Q{H + Q)) - Q{H)]bb\ 



(2.101) 



(2.102) 



(2.103) 



(2.104) 



unless n{H + n) = n{H). 

In general, it is not possible to solve equations ( |2.93D except for some trivial models. One 
can, however, obtain approximate solutions by means of perturbation theory. If we write H = 
Hq + XH' and expand 



3=0 j=0 



(2.105) 
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then we obtain the coefficients Vlj and bj from [|33|J3^] 



[Ho, ilj] 

[Ho, bj] + VLobj 



H'], 

[bj-i,H']-j:i=i ^kbj-k 



(2.106) 



For concreteness and simplicity we consider dimensionless anharmonic oscillators ( |1.29D . 
In two recent papers Fernandez |3^ ] obtained Qi(Ho) for several values of m, and Qi(Ho) 



and Q2{a, a^) for m = 4. These operators enable us to derive energy differences in the form of 
A-power series: 



= E ^n,jX', En = J2 En,jX' (2.107) 
j=0 j=0 

In the present case fio = 1 so that fii commutes with Ho- For sufficiently small A we have 
b = a + 0{\) and we may choose m = n + 1 in equation ( [2.96[ ), so that 



Q ■ = E ■ — E 



For the first order we have: 



1 =< n\Qi\n >, 



(2.108) 



(2.109) 



where |?7.) is an unperturbed state. The operator is diagonal and does not contribute to the 
correction of second order; therefore 



^n,2 =< n\Q2\n >; . 



(2.110) 



All the expressions of first and second order derived by means of the method just outlined \ p5 



already satisfy equations ( [2.108[ - P7l 10] ). Fernandez has obtained fii for different m in this 
approach. Here we establish a general form of ^7l valid for all integer values of m. To do so we 
use the fact that ^7l is always a function of the unperturbed Hamiltonian Ho- Now 



< n|f2i|?T, > 



En+l — En 

.(f-1) 



1 + 

T rrYi ^ O O ^ 



m— 2r ^ 



1 v-(f-l) 

— T7T 



E\ 2 '4- m/^ Ta—lrr^ 
r=0 '^2r ^2r ^. 



2 

^ ^ I 1 Y^^"2"~^^ J. mr^ m~2r 



■(n- 



.iri^+l)!(2L|2r__l)! 



' m—2r I 



r(ata+l) 



' 2^m 



■ r(ata 



(f-l) 



m — 2 r 

2 . 



' m— 2r ^ 



n > 



L2r ^2r O 77t— 2- 



r(Ho 



+2)! ^ 2 



Therefore we have 



i'l — "TH / , r2r ^2r O m-2r 

2~m 



(f-l) 



r(/7o + ^) ,m-2r. 



r(i7o-^ + i)! 



(2.111) 



(2.112) 



r=0 V-'-'U 2 '2^ 

This coincides exactly with the first order frequency operator obtained by other techniques. 
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2.8 Comparison among the different approaches 



We have seen that different approaches may be used to obtain the first order frequency operator 
and the first order position operator. But all the results obtained in these methods are not the 
same. In this section we have checked the equivalence of different methods and have established 
some relationships between apparently different but equivalent frequency operators. To begin 
with, let us compare the first order cases in the next subsection and higher order cases in the 
subsequent subsection. 



2.8.1 First order case 



In the preceding sections we have obtained the same frequency operator from perturbation 
theory in the interaction picture, renormalization group, near-identity transformation, and eigen 
operator techniques. In these techniques we usually write 



X{t) = A exp{-int)a{0) + a\0) exp(ifit) 



(2.113) 



where Vl = ^1{Hq) = 1 + \^Ii{Hq) + 0{X^) and A is the normalization factor. However, this 
form of the frequency operator appears to be different from the one derived from the Taylor 
series method [|3^, |3^] and multiple- scale analysis [ pO| , pT] and |3^]. In these approaches position 
operator X(t) is usually expressed as 

X{t) = cos(tut)X(0) + X(0) cos(cjt) + sm{ujt)X{0) + X{0) sin(wt) + 0{\), (2.1 14) 

where to = uj{Hq) = 1 + \uji{Hq) + 0{X^). Taking into account a particular form of equation 
( ^JWj ): aF{Ho) = F{Hq + l)a, F{Ho)a^ = a^F{Ho + 1), we easily rewrite equation (^JT§ 
as 



X{t) 



cos At 



Ui — tUi 



exp 



exp it + iXt 



-it — iXt 
+ 0(A2) 



(2.115) 



where iu[ = lui (Hq + 1) . Upon comparing equations ( |2.113| ) and ( |2.115| ), and taking into account 



that cos 



(At- 



1 + O(A^), we conclude that 



(2.116) 



The inverse transformation 



Co — 1; Cj 



2 



(2.117) 



also gives a suitable analytical expression because both Vli and ui are polynomial functions of 
Ho. 
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2.8.2 Corrections of higher order 



The calculation of perturbation corrections of higher order and comparison of expressions com- 
ing from different methods is considerably more difficult. We can, however, draw some useful 



conclusions from available results. For example, Aks and Carhart p9| ] showed that the Hamil- 
tonian operator H(p{t), x(t)), in terms of their perturbation expressions free from secular terms 
for x{t) and p{t), is diagonal and gives the correct energy through first order perturbation theory. 
Unfortunately, they did not verify whether their second-order results were consistent with this 
criterion. 

Egusquiza and Valle Basagoiti obtained an expression of the form 



a{T] = aexp<; — + 
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exp 



-3trA(Ho + |) 



+ 



5 H„- 



-1 



(2.118) 
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by means of renormalization [ pBQ . The second-order part of their frequency operator 



64 V'''^ + 



yields 

which does not agree with the well-known result 



n - linf^ln - 1) = - — (5^2 - 1) 
12 1 / g^v ; 



En,2 — E,, 



^^-l,2 
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(2.119) 
(2.120) 

(2.121) 



that one easily obtains by means of Rayleigh-Schrodinger perturbation theory p5|] 



By means of the method outlined in Section 2.8 Fernandez obtained [|36|] 



^2 



_3_ 

16 V 4 



-^(69i72 + 5mo-12X^ + if) 



(2.122) 



which is consistent with equations ( [2.110[ ) and ( |2.121D . This expression shows off-diagonal terms 
and agrees with the perturbation expansion of Speliotopoulos's frequency operator equation pll ] 
when A = 4e. There seems to be a conceptual oversight in Egusquiza and Valle Basagoiti's 
correction of second order to the frequency operator because it does not show offdiagonal terms 
and does not give the correct energy difference P3|]. To be more precise, the second order 
frequency operator obtained by them differs from all other existing results [see equation 2.10 of 
ref. ^ and the discrepancy is not due to a difference in operator ordering. This fact may be 
established from their expression for the energy difference: 



^ / 3A 3A^ 2 



1^ 



0(A3 



(2.123) 
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They obtained ( [2 .1231 ) as a consequence of equation 2.10 of ref. Equation ( |2.123| ) is equa 



tion 2.11 of ref. p3| ] and it differs from the energy difference obtained by simple perturbation 
calculation which is 



En — En- 
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(51n2 + 21) + 0(A' 



(2.124) 



Possibly the relevant secular term chosen by them (see equation 2.8 of []33|]) is incomplete and 
this causes the discrepancy in the expression for frequency operator obtained by them and others. 
By means of a completely different approach we obtain 



UJ2 



which is consistent with their recipe for the second-order energy difference [ p7| ] 



2 {En,2 - En^i,2) = {n\uj2\n) + {n- l\uj2\n - 1). 



(2.125) 



(2.126) 



Equations ( |2.121| ) and ( |2.126| ) suggest that the diagonal parts of the operators ( |2.122D and ( |2.125D 
are related by 



(2.127) 



as one easily verifies by means of the operators themselves. 

Aks and Carhart [ p9| ] forced the frequency operator to be diagonal at all perturbation orders. 
They proved that their results give the Rayleigh-Schrodinger perturbation series of first order 
where we know that the frequency operator is diagonal and correctly given by all the approaches. 
We have seen that from second order onwards the frequency operator is a polynomial function of 
H, and, therefore, it should exhibit off-diagonal terms for the unperturbed states. It is not clear 



to us if the results of Aks and Carhart [ P9| ] are correct at second order because they did not give 
the corresponding expression of the coordinate explicitly. Moreover, it is difficult to compare 
the results of Aks and Carhart with others because the way they split the frequency operator into 
the different corrections to the coordinate operator differs markedly from the other techniques 
discussed here. 

Here we can also check that if we put H = Hq + ^x^ in the expression of \{H) obtained by 
Speliotopoloulos pT] ] and take terms up to the second order in A (where A = ^) then we obtain 
the second order frequency operator obtained by Fernandez [|3^]. So the results of Fernandez 
and Speliotopoloulos are in exact accordance with each other. Actually according to the work 
of Speliotopoloulos frequency operator is a function of the total Hamiltonian H = Hq + XV. 
Now if we express the frequency operator in power series of A as = i^o + Ai7i(ifo + ^V) + 
\^fl2iHo + W) + .... then we observe that the first order correction to the frequency operator 
(Vli) should be a function of Hq only since even the lowest power of V will give terms of the 
order of A^ . So it is expected that the second order frequency operator will be a function of 
Hq and V but no cross term will be present but third and higher order frequency operator will 
contain cross terms like H^V. This is really in accordance with the available results and we can 
recall that the second order frequency operator for the quartic oscillator obtained by us is 
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{69H^ - I2x^ + 
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This matches with the argument given above. 

Thus we conclude that apparent differences among perturbation corrections of first and sec- 
ond orders to the frequency operator derived by different alternative approaches are simply due to 
different arrangements and ordering of products of noncommuting quantum-mechanical observ- 
able. In some cases it is not difficult to obtain the relationships between the different forms of the 
frequency operator as shown above. Therefore we can use any of these solutions to investigate 
the quantum fluctuations of coherent light in nonlinear media. 
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Chapter 3 

Phase fluctuations of coherent light coupled 
to a nonlinear medium of inversion 
symmetry 



The quantum phase problem may be stated as: How can one write down a quantum mechanical 
operator corresponding to the phase of a harmonic oscillator or equivalently a single mode of 
electromagnetic field? The concept of phase plays a very important role in the understanding 
of basic physics. So people are trying to answer this question from different points of view and 
are debating on this issue since Dirac p2| ] initiated the search of a quantum mechanical operator 
corresponding to the phase of a harmonic oscillator in 1927. Due to this seventy-five years long 
search and healthy debate an extensive amount of literature on this topic is now available but a 
satisfactory Hermitian phase operator is still unavailable. However, the observable nature of the 
phase demands that the corresponding operator should be a Hermitian one. There are many good 



reviews on this topic 058| - |62| ]. Some of them are written in a lucid manner [ p8l , pl| ] and some are 
written from much more mathematical and technical points of view [|62[]. 

A complete description of a harmonic oscillator is possible from the idea of a phase. Again 
the interference and diffraction pattern of electromagnetic waves can be explained convincingly 
with the help of a phase concept. Thus the phase in classical physics is well understood for a 
long time. We do not have such resemblance of phases in quantum mechanics. But it is widely 
believed that the quantum phase is one of the distinguishing features between quantum and clas- 
sical physics. Actually quantum phase is responsible for all quantum mechanical interference 
phenomena. Even the discrete eigenvalues in quantum theory can be viewed as a quantum phase 
condition for the Schrodinger equation [^. The examples of quantum phase are from double 
slit experiment [ |53| ] to the dark resonances and from Aharonov-Bohm [ |55| ] effect to Bose- 
Einstein condensate [|6^]. 

Quantum mechanical phases can be classified into two main parts depending upon their ori- 
gin. The phase having geometrical origin is quite familiar to us and is called geometric quantum 
phase (e.g Berry's phase [|57|], Aharonov Anandan phase The other type of phase arises 

due to the dynamical nature of the system. Geometric phase is discussed in chapter 6. In this 
chapter we briefly review the quantum phase problem and study the quantum phase fluctuations 
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of coherent light coupled to a third order nonlinear medium of inversion symmetry. 



3.1 The quantum phase: Dirac approach 



Dirac started his work with the assumption that the annihilation operator a can be factored out 



into a Hermitian function f{N) of the number operator N and a unitary operator U [^. The 
later one (i.e.U) defines the Hermitian phase operator (p in the following way 



U = expiicp). 



Hence in the Dirac formalism 



a = exp{i(j)) f {N) . 
Now the functional form of f{N) may easily be obtained by using 

N = a)a = f. 
Thus the explicit expression for a is given by 

a = exp{i(t>)N^ 



(3.1) 
(3.2) 

(3.3) 
(3.4) 



where equations ( p.2D and ( p.3[ ) have been used. The expression ( P^ looks like the polar de- 
composition of a complex number (a). The usual commutation relation between a and is 

aa^ — a) a = 1. 



(3.5) 

(3.6) 
(3.7) 

(3.8) 

The justification of the above condition ( |3.8| ) may be simply given in the following way. If 
equation ( p.Sj ) is true then from the method of induction we have 



From equations ( p4| ) and ( |3.5D we have 

exp{i(f)) Nexp{—i(j)) — N = 1 

or, 

exp{i(f))N — Nexp{i(f)) = exp{i4 



Equation (|3.7[) is valid only if the following commutation relation holds good 



[N,4>]=z. 



[N, 0"] = 



n-l 



(3.9) 



Multiplying both sides of equation ( |3.9D by ^ and summing from to oo we arrive at equation 
(^77|). Now the uncertainty relation corresponding to equation ( |3.8D is given by 



AN Ad) > 



(3.10) 



The equations ( |3.8| ) and ( |3.10| ) have many problems. Let us see them one by one. 
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3.1.1 Problems with the Dirac approach 

Problem- 1: The uncertainty relation ( |3.10| ) allows the uncertainty in </> (i.e A0) greater than 27r 
for AA^ < This is impossible and is inconsistent with the sense of a phase. This problem is 
known as the periodicity problem. In this context we can recall that the uncertainty relation for 
Lz and (p also has the same problem [^^-[7^ . 

Problem- 2: The Hermitian nature of (as it is expected in quantum mechanics) demands that 
the operator U should be an unitary one. However U is not an unitary operator. 
Problem 3: In 1963 Louisell [ [7T| ] showed that equation ( |J!8| ) leads to confusing results when one 
attempts to evaluate the matrix element of in a representation in which is diagonal. Since, 

< l\N(t) - (f)N\n >= i <l\n> 

so we have 

<mn>=-tj^^ (3.11) 
[n-l) 

where N\n >= n\n > has been used. If the eigenvalues of N are integer then the right hand side 
of equation ( p.l ip becomes undefined. Only in the limit of large values of n and / (Correspon- 
dence limit) 4> takes on a definite meaning. Thus has a definite meaning only in the classical 
limit. This is quite consistent with the classical idea of phases. 



3.1.2 Why does the Hermitian phase operator not exist? 

In the previous section we have seen that the Hermitian phase operator does not exist. But a 
question remained unanswered: Why does a Hermitian phase operator not exist? J. Sarfatt, one 
of the pioneer in quantum phase problem, addressed this question for the first time. He wrote 



a nice paper in connection with the complexity of the introduction of quantum phase [^. In 



Sarfatt words [|72[]'... Recent work by L. Susskind, J. Glogower and J. Sarafatt shows that it is 
impossible to define a phase operator because of the existence of a lowest state for the number 
operator of the oscillator. Thus the uncertainty relation ANAcp > 1 is meaningless....' 

The main cause for the nonunitarity of U is that the number state matrix representation of the 
diagonal number operator is bound from below. So a bare phase operator does not exist. This 
observation provoked people to introduce geometrical functions of phase operators. 



3.2 Periodic function of phase: Louisell 's approach 



Louisell [ j7T| ] proposed a self consistent way to avoid these difficulties. Actually we can write 
equation ( p.9D in an alternative form as 

[N,r]=^nr^' = ^^. (3.12) 
Therefore for any polynomial function of Lp we have a commutation relation 
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[N,P 



.dP{(p) 



(3.13) 



Louisell started his work with the fundamental commutation relation ( [3.13D and supposing that 
P(0) is a periodic function of period 27r. He considered sin(0) and cos(0) as the Hermitian 
operators which satisfy the following commutation relations, 



[iV, cos(0)] = —ism{(f)) 



and 



[A^, cos(0)] = i sin(0). 
Therefore, the uncertainty relations take the following forms 

AiVAcos(0) > ^1 < sin(0) > 
AA^Asin(0) > ^1 < cos(0) > 



(3.14) 
(3.15) 

(3.16) 
(3.17) 



It is to be noted that the above uncertainty relations ( |3.16| and [3.17D are free from problem 1 . 
Instead of a bare phase operator, Louisell defined a meaningful phase operator in terms of the 
periodic function of phase. But the following questions remain unanswered: 

1) What are the explicit forms of sine and cosine operators? 

2) Are they Hermitian? 



3.3 sine and cosine operators: Susskind and Glogower approach 



The questions asked in the previous subsection were properly addressed by Susskind and Glo- 
gower (SG) [[73|]. They reaffirmed that London's result JT^ , [751 ] (i.e U is not unitary) are indeed 
correct and obtained an explicit form of sine operator S and cosine operator C as 



and 



S 



c 



-(E- 



1 

2i 



{N + l)-2 



-a — a' 



(AT + 1)2 



2^ 



+ ^4 



-a + a 



(3.18) 



(3.19) 



XN+1)2 {N+l)2 

It is easy to check that the operators S ( |3.18| ) and C ( p.l9D are Hermitian. Therefore, SG approach 
is fully capable of giving the answers to the questions remained unsolved after Louisell's work. 
Now the commutator of C and S may be calculated by using (p.l9| and p.l8[) and is given by 



[C, S] = -P' 



(3.20) 



where P° is the projection onto the ground state. The noncommuting nature of C and S ( |3.20| ) is 
one of the serious drawbacks of SG approach. This is because the cosine and sine functions of a 
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real phase should commute. Moreover, the sum of the squares of the S and C operators is given 
by 

C^ + S^ = 1- ^P°. (3.21) 

Therefore, in SG approach we got satisfactory answers to the previous questions. However, 
we ended up with two new problems which essentially go against our perception. 



3.4 Measured phase operators: Pegg and Barnett approach 

In 1986 Pegg and Barnett [|7^] proposed a Hermitian phase operator in an extended Hilbert space 
containing negative number states. Such a proposal seems to be unphysical. So very soon they 
gave up this proposal and started their new venture with a new proposal. They proposed that the 
infinite dimensional harmonic oscillator Hilbert space should be truncated to (s + 1) dimensions 
and expectation values of the quantum mechanical operators should be taken in this truncated 
Hilbert space only. After doing that one should take the limit s — > oo. This prescription of 
working in a truncated Hilbert space was not new in Physics. Popov-Yarunin, Sinha and many 
others already did some work on the finite dimensional Hilbert space. But the Pegg-Bamett 
prescription made the problem much more user friendly and it also defined a Hermitian phase 
operator in the truncated Hilbert space. So this formalism got a lot of attention of the people. 
Pegg and Barnett started their work by defining a state having a well defined phase p3| - p4| ] 

as 

s 

\e >= lims^oois + 1)"^ exp(m6')|n > (3.22) 

n=0 

where \n > is the number state spanning the (s + 1) -dimensional space. In this space they 
selected a reference phase state 

s 

1^0 >= (s + l)-5 ^exp(m^o)|^ > (3.23) 

n=0 

and found the subset of states \6m >, defined by replacing Oq by 9m in equation ( [3.23[ ), which are 
orthogonal to this state. Thus we have 

s 

< OmlOo >= (s + l)-i exp(m(eo - OJ). (3.24) 

n=0 

One can easily prove that the phase states \9o > and \9m > are orthogonal to each other if the 
following relation holds 

9m = eo + —-{m = 0,l,....s). (3.25) 
s + 1 



The state (3.23) is an over-complete one and \9m > forms a complete set of orthogonal basis 



vectors spanning the state space. Therefore, the number state \n > can be expanded in terms of 
1 6*^ > as follows 

s s 

\n >= Y \^rn >< 9m\n >= (s + 1)"^ Y exp{-in9m)\9m > ■ (3.26) 

m=0 m=0 
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In this approach phase and number basis states are so related that a system in a number state is 
equally likely to be in any state \9m > and a system prepared in a phase state is equally likely to 
be found in any number state [0]. Now the problem is to find a unitary operator expg{i(f)) whose 
eigenstates are the phase states \9rn >'■ 



eycpg{i(f))\9m >= exp{i9m)\0m > 



(3.27) 



>= exp{-i9,rn)\0m > ■ 



(3.28) 



Here, the caret indicates that the whole expression is an operator. From equation ( |3.26| and |3.27D 
we have 



ex^g{i(f))\n >= [s + 1)-^ expH(n - l)9„Mm >= \n - 1 > . (3.29) 

m=0 

For the vacuum state (n = 0) the resulting state is 

(s + 1)-^ J2,^exp{i9^)\9^ >= (s + l)-^exp[i(s + l)6'o] Em exp(-is6l„) |6'^ > 



exp[i{s + l)9o]\s > 



. (3.30) 



Now from equations ( p.29| and |3.30| ) one can see that the number state representation of expg(i^ 
is 



expg[i^ 



|0 >< 1| + |1 >< 2| + + |s - 1 >< s| + exp[i{s + l)9o]\s >< 0| (3.31) 



and expg(—i(f)) is just the Hermitian conjugate to it. Thus we have unitary operator expg{i(f)) and 
we can define cosine and sine operators in terms of that. Those operators have properties that 
are much more desirable and familiar for the description of phase. In particular in Pegg-Bamett 
(PB) formalism we find 

[cose0]2 + [shTe^]' = 1 (3.32) 



[cos6)0, smQ(f)] = 



(3.33) 



where 



and 



< n|[cos6)0]^|n >=< n| [singf^j^ln >= ^'^^^^ 



= -[expg{i(j)) + expg{-i(f 



sine0 = —[expg{i(f)) - expg(-i0)]. 



(3.34) 
(3.35) 
(3.36) 



One of the main differences between SG formalism and PB formalism is that the action of the 
exp{i(f)) operator on the vacuum gives different results. This is basically responsible for all the 
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differences in operator relations we have obtained in two formalisms. From equation ( p.34D it 
is clear that in PB formalism the definition of phase operator is consistent with the fact that the 
phase of the vacuum is random but the case is not so in the S -G formalism. 

Though in Pegg-Bamett formalism (p has the desirable properties of phase but still this ap- 
proach is not free from objections. In this approach one measures the physical quantities at first 
in an arbitrarily large but finite state space of s + 1 dimension and then takes the limit as s tends 
to oo. These two operation do not commute so their applicability is questionable []77|]. But Pegg 
and Barnett demand that their approach is physically indistinguishable with the conventional in- 
finite state space model p3|]. Correctness of the results obtained in different approaches may be 
compared only by experiments. Interesting experiments on quantum phase have started coming 
in@. 

People, however, use both Susskind-Glogower (SO) and Pegg- Barnett (PB) [§, ^ 

formalisms for the studies of phase properties and the phase fluctuations of various physical sys- 
tems. For example, the phase fluctuations of coherent light interacting with a nonlinear nonab- 
sorbing medium of inversion symmetry are reported in the recent past [|^J^]. It is found that both 
SG and PB [||] formalisms lead to same type (qualitatively) of phase fluctuations. Keeping 
these facts in mind, we have chosen PB approach for the present investigation. 

In chapter 1 we have seen that the interaction of a single mode of electromagnetic field with 
a third order nonlinear medium can be described by the Hamiltonian ( |1.31[ ). The equation of mo- 
tion corresponding to the Hamiltonian ( |1.31D is given by equation ( |1.32D which contains cubic 
nonlinearity in field operator X. The exact analytical solution of the equation ( |1.32| ) is not avail- 
able. However, under RWA the model Hamiltonian reduces to an exactly solvable form ( |1.33D . 
Gerry and Lynch |Q] used this RW approximated Hamiltonian ( |1.33| ) to study the quantum 
phase fluctuations. On the other hand, different techniques to obtain approximate solution of 
( I1.32D is discussed in chapter2. Here we use first order operator solution of ( |1.32| ) obtained 
in chapter 2 to study the quantum phase fluctuations of coherent light in third order nonlinear 
medium. 



3.5 Time evolution of the useful operators 

The purpose of the present section is to calculate the phase fluctuations of coherent light interact- 
ing with the nonlinear medium of inversion symmetry. The corresponding Hamiltonian is given 
by equation ( |1.31| ) and the fluctuations are measured with respect to the initial coherent state 
I a > which is defined as the right eigenket of the annihilation operator a(0) corresponding to the 
eigenvalue equation a(0) |a >= a\a > [^^. The eigenvalue a is in general complex and may be 
written as a = \a\exp{i9), where 9 is the phase angle of a. The modulus square of a {i.e\a\'^) 
gives the number of photons present in the field prior to the interaction. From chapter 2 we can 
write the time evolution of the annihilation operator a{t) as 

a{t) = X{t) + iX{t) 

= Dia(O) + Dsa^O) (3.37) 
- Dsa^iO) + D^a^^O) + D5a^\0)a{0) + DGa^{0)a\0) 
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where the parameters 





= (1 


D2 


= -i- 


D3 


= 






D, 


_ ,-3A 


and 






= 



:3A. 



3A 
4 



t) exp(— it) 



ij sint exp(— 2it) 
i| sin 2texp(zt) 
sint 



(3.38) 



The time evolution of the creation operator a^{t) is the Hermitian conjugate of ( 13.37D . Thus the 
number operator is 

N{t) = a^{t)a{t) -- 



at(0)a(0) + (D*D2at2(0) + h.c) 

- (DlD3a^{0)a^{0) + h.c) - (DlDia^^{0) + h.c) (3.39) 

- (DtD5a"^^(0)a(0) + h.c) - (D*D6a"^^(0)a2(0) + h.c) , 

where h.c stands for the Hermitian conjugate. The parameters D^, D^, D^, Dl, and are 



the complex conjugate of Di, D2, -D3, -D4, -D5 and Dq respectively. The equation ([3.39| ) is used 
to obtain 



N\t) = (at2(0)a2(0) + at(0)a(0)) 

- 2DID3 (at2(0)a^(0) + 2at(0)a3(0)) + h.c 

- '2DID^ fat5(0)a(0) + 2a^^{0)) + h.c 

- '2DID5 (a^\0)a\0) + 2at3(0)a(0)) + h.c 

- '2DID6 (at3(0)a3(0) + 2at2(0)a2(0)) + h.c 



2DID2 at3(0)a(0) + at2(0) + h.c 



(3.40) 



The terms beyond the linear power of A are neglected in equation ( [3. 40] ). 

The exponential of phase operator E and its Hermitian conjugate E^^ under the PB formalism 
are given by 

E = fiV+i' 



2 1 '^(^) 



(3.41) 



E^ = {N+iy 

where is the average number of photons present in the radiation field after interaction. Using 



( 3.39|) N can be written as 



N 



\a\ 



1 + ^ (2(3 + 2|a|2) sintsin(t-2^) + |a|2sin2tsin2(t - 20) 



(3.42) 



Interestingly, N depends on the coupling constant A, phase angle 9 and on the free evolution 
time t. Thus the number of photons are not conserved. The result is not surprising since the 
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nonconserving energy terms are included in the model Hamiltonian ( |1.31D . However, in earlier 
studies the photon numbers were conserved. 

The usual cosine and sine of the phase operator are defined in the following way 



C 

s 



(3.43) 



The expectation values of the operators C and S are given by 



<C> 



^'N+\] '''[{Di + Dl)a + {Dl + D2)a* -{D^ + DDa" 



~{D* + Di)a*'^ - (D5 + Dl)\a\'^a* - [D^ + Dq) \a\'^a 

2 



( iV + ^ ) " [{Di - Dl)a - {Dl - D2)a* - {D3 - Dl)a^ 



(3.44) 



<S> ■■ 

+ {Dl - Di)a*'^ - (L>5 - Dl)\a\'^a* + (D* - Dq) \a\'^a] 
where the equations ( [3.37| - p39| ), and the equations ( |3.41[ - |3.43| ) are used. Again, the square of the 



averages are 

<C>^ 



-1 



l{N + l) [{{Df + 2DiD*)«2 ^ ^^^1 ^ 2 + (D,D2 + c.c)} |«| 



-2 ({(Z^iZ^3 + D^Dl)a^ + c.c} + {{D,D* + D^D^)\a\^a*^ + c.c} 
+ + D^Dl) + c.c} \a\^ + {{D^Dl + DiD^)\a\''a^ + c.c})] 



(3.45) 



-\ (n + I) [{{Dl - 2DiD*)«2 + ^^^1 _ 2 _ {^D.D^ + c.c)} |a| 

-2 ({(D1D3 - + c.c} - {{D^Dl - D^D^)\a\^a*^ + c.c} 

V {{D,D, - D,D;) + c.c} \a\^ - {{D.D; - DM\a\^a^ + c.c})] 



(3.46) 

where c.c stands for the complex conjugate. Using equations ( |3.39| - |3.46| ) the second order vari- 
ances of C, S and can be written as 

(AC)2 = l(N + iy\{\Di\^ + (DiD2 + c.c)}-3{{DiD* + DiD4)a*^ + c.c} 
- 2 + D,D*,)\a\^ + c.c)} - {{DiDl + + c.c}] , 

(3.47) 



-\ + I) ' K-I^il' + (^1^2 + c.c)} + 3 {{D^D* - D^D^)a*^ + c.c} 
2 - D^DDlal-" + c.c)} + {(Z^iD* - D,D,)a^ + c.c}] , 



(3.48) 



(AiV)^ 



(1 + A [(3 + 4|an sin t sin(t - 2^) + |ap sin 2t sin 2(t - 2^)]) 



(3.49) 
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3.6 Phase fluctuations 



The usual parameters for the purpose of calculation of the phase fluctuations (for a fixed value of 
A) are defined as [BIS] 



U (e, t, \af) = (AiV)2 [{AS)^ + (AC) 



<S>' + <C>' 



(3.50) 



and 



(3.51) 



Q (e,t, = S {e,t, /< C 



(3.52) 



Now we can, in principle, analytically calculate f/, S and Q by using equations ( [3.42| - p.49| ). The 
equations ( |3.50| - [3.52| ) assume the following forms. 



u {e,t, \a\ 



1 + f (6(1 + 2|a|2) sintsin(t -6) + 3|a|2sin2tsin2(t - 26)) 



(3.53) 



^1 



h t, |a| 



= ^(iV+i)"'[l + f {6(l + 2|«|2)sin2t 

+ 6\a\Hsm2{t- 6) +3\ sin 29 sin 2t 

+ 4(3 + 2|ansintsin(t-2^) + 4|a|^sin2t sin2(t - 2^)}] 



(3.54) 



and 



Q{9,t,\a\ 



+ 



4cos2(t_e) 



1 + f {6(1 + 2|a|2) sin^ t + 4(3 + 4:\af) sintsin(t - 29) 



\a\ 



■ [6t sm2{t-29) + 3 sin 2^ sin 2t + 4 sin 2t sin 2 (t - 2^)] } 



- scos^Ve) {-^tsm2{t- 9) + (6 + 4|an sint sin(t - 2^) 

- 6(1 + laHsin^t - 2|ap sint sin(3t - 49) + |apsin2tsin2(t - 29) 

- lap sin 2t sin 2^ - 6|aptsin2(t - 9)}] . 



(3.55) 

Hence the equations ( |3.53| - p.55| ) are our desired results. In the derivation of the equation ( p.55p , 

-1 



(I 



a\ 



and Qc 



are the initial 



we assume 7^ 0. Now, f/o - 2, - 4i"i (^i-i . 2j — - 4cos2(t-0) 
(i.e A = 0) values of U, S and Q respectively. Thus Uq, Sq and Qq signify the information 
about the phase of the input coherent light. The suitable choice of t may cause the enhancement 
and reduction of all the above parameters compared to their initial values. It is to be noted that 
the parameters S and Q contain the secular terms proportional to t. However, it is not a serious 
problem since the product At is small [|32|]. The equations ( |3.53| - p.55| ) are good enough to have 
the flavor of analytical results. Now we would like to discuss two interesting special cases. 
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3.6.1 The vacuum field effect 



The radiation field with zero photon is termed as the vacuum field. The interaction of vacuum 
field with the medium gives rise to some interesting quantum electrodynamic effects []5^]. The 
purpose of this subsection is to study the effects of the vacuum field on the phase fluctuations of 
input coherent light. In case of vacuum field = 0) the equations (p.53|) and (p.55|) reduce to 



U{9, t) = f/o f 1 + y sin t sin(t - 9) 



(3.56) 



3A 



Q{e,t) = Qo [l + f sint{sint + 2sin(t-2g)}- 

X {-tsin2(t -9) + smtsm{t - 29) - sin^t} 



(3.57) 



where |ap is used to derive ( |3.57| ). Interestingly, the vacuum field itself couples with the 



medium and gives rise to the condition A 7^ 0. Thus the phase fluctuations for vacuum field are of 
purely quantum electrodynamic in nature. Now for = 0, the parameters U and Q are enhanced 
compared to Uq and Qo respectively. The corresponding maximum fluctuation of U is obtained 
if t is an odd multiple of n/2. The value of Q is infinite if t becomes an odd multiple of n/2. 
For 9^0, the parameters U and Q may be reduced or enhanced by the suitable choices of t. It 
is to be noted that the parameter U is 0.5 and is independent of 9 in the earlier occasions [^p. 
In case of the present work, however, the parameter U depends on 9 and t. The parameter S is 
identically zero and coincides exactly with the earlier results 

3.6.2 Phase of the input coherent light = j 



The equations ( p.53[ - P35| ) are still complicated even for 6* = 7r/4. A further simplification is 



made with the choice t = 7r/4. The equation ( [3.53D reduces to a simple form 

A huge reduction of U is possible with the increase of the photon number \a\'^. However, care 
should be taken about the condition of the solution during such increase. The circular nature 
of the trigonometric function ensures the occurrence of similar reductions for other values of t 
(t = 2mn + I where m is an integer). Now for t = n/2, we have 

C/(f.^.H')=C„(n-5^{l + 2Hn). (3.59) 



The equation ( [3.58D clearly shows the enhancement of U parameter as the photon number 



a 



2 



increases. Thus we conclude by noting that the parameter U may decrease ( p.58D or increase 
( [3.59D with the increase of |ap by suitable choices of t. Similarly, one can obtain the reduction 
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and enhancement of the remaining two parameters (Q and S) by suitable manipulations of free 
evolution time. 

In the earlier works [^p, the parameter Q is found to decrease with the increase of pho- 
ton number till a minimum is reached. Subsequent increase of lap causes the increase of Q. 
However, the parameters U and S are always enhanced compared to their initial values as |ap 
increases. Hence, the present results are in sharp contrast with those of the earlier studies 
Hi . The above differences are attributed due to the fact that the nonconserving energy terms are 
included in the model Hamiltonian. 



3.7 Conclusion 

The phase fluctuation of coherent light interacting with a nonlinear medium of inversion sym- 
metry is carried out by using the PB formalism. The usual parameters for this purpose are U, Q 



and 5". These parameters ( p.53hp.55| ) are found to depend on 6*, \a\ and t. It is interesting to note 
that the free evolution time t is absent in the earlier works [^p. However, the presence of t is 
automatic in the present calculation. It accounts for the fact that the nonconserving energy terms 
are present in our model Hamiltonian. 

The effect of vacuum field on the phase fluctuation parameters are expressed in closed ana- 
lytical forms. It is found that the enhancement of U and Q are possible when the phase angle 
= 0. However, for 9^0, both reduction and enhancement of U and Q are possible by suitable 
choices of t. The observed results are in sharp contrast with those of the earlier studies [^H]. The 
parameter S reduces to zero for vacuum field and agrees exactly with those of the earlier works 



Apart from the general expressions for U, Q and S, we also made a qualitative comparison 
between the present results and the results obtained earlier for the identical physical system. 
For 6 = IT /A, we obtain the reduction and enhancement of phase parameters (U, Q and S) by 
suitable manipulation of free evolution time t. Those results are in sharp contrast with the results 
already obtained for the same physical system [^§p. Clearly, the free evolution time t makes the 
differences. 



Of late the preparation of quantum states have been reported by several laboratories [ pl| - |830 . 
These production of quantum states have opened up the possibilities of experimental studies on 
quantum phase and hence the verification of the present results. 
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Chapter 4 



Squeezing of coherent light coupled to a 
third order nonlinear medium 



According to the Heisenberg's uncertainty principle we can't measure position (X) and momen- 
tum (X) of a particle with greater accuracy than 

AX AX > ^ (4.1) 

where uncertainty (^C*) in 

any arbitrary variable C is defined as (AC) = {< > — < C >^)^. Although, there is no 
position operator for photon, fundamental canonical variables for light (component of electric 
field and magnetic field) satisfy an uncertainty relation of the form ( [4.ip and an electromagnetic 
field is said to be electrically squeezed field if AX is less than the vacuum field uncertainties 
in the quadrature (i.e. (AX)^ < |). Correspondingly a magnetically squeezed field is one for 
which (AX)^ < i. Here we are not going to discuss the phenomenon of squeezing in detail. So 



for a review on squeezed states, one can see the references [ p4| , p5[ |. In this chapter we want to 
examine the possibilities of observing squeezing in a physical situation in which a single mode 
of the quantized electromagnetic field, initially prepared in the coherent state (coherent light), is 
allowed to interact with a nonlinear nonabsorbing medium of inversion symmetry. 

4.1 Application of quantum quartic oscillator: The squeezed states. 

We have already shown that for the lowest order of nonlinearity (i.e cubic nonlinearity), corre- 
sponding Hamiltonian of the system is given by the Hamiltonian of a quartic anharmonic oscil- 
lator ( |1.31D . Equation of motion of the corresponding oscillator is ( |1.32D and the second order 
operator solution of ( 11.32| ) is simply given by ( |2.25[ ). 



Let us now express the quadrature operator X obtained in ( |2.25| ) in terms of the initial {i.et = 
0) annihilation and creation operators. It is given by 

X(t) = [Eia(O) +^20^(0) + E3at2(0)a(0) + ^40^(0) 

(4.2) 



E5at(0)a4(0) + E6at2(0)a3(0) 



+ h.c 
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where h.c stands for the Hermitian conjugate. The remaining quadrature operator X{t) may 
simply be obtained by taking the time derivative of (|4.2|). The parameters Ei [i = 1, 2, 3, 4, 5 and 
6) are complex and are given by 



El 



+ 
+ 



cost — isint — ^ {tsint + i{tcost — sint)} 
^ { (468t sin t - 63 cos t + 63 cos 3t - 2 1 6t^ cos t) 



1 

512 



;il88t COS t - 1053 sin t - 45 sin 3t + 216^^ sin t) }] 



(4.3) 



E, = ^ 



1 



-j^ {(cost — cos3t) + i(sin3t — 3sint)} 



+ 
+ 



^ {(156 cost - 192t sin t - 156 cos 3t - 144t sin 3t) 



i(156 sin 3t - 324 sin t - 156t cos 3t) }] 



E, = 4. 



1 

V2 



+ 
+ 

E. 



^ { (936t sin t - 126 cos t + 126 cos 3t - 432t2 cos t) 



— ^ {tsint — i(tcost — sint)} 



and 



z(2376t cos t - 2106 sin t - 90 sin 3t + 432t2 sin t) }] 

= 256^2 ['-^ ^ ~ ^'^^ sin t — 6 cos 3t + cos 5t) 
+ z(6 sin 3t — sint — 12t cost — sin 5t)] 

E^ = ^^[(39cost-48tsint-39cos3t-36tsin3t) 
+ 2(39 sin3t — 81 sin t — 36t cos3t)] 

Eq = ^g^[(156tsint-21cost + 21cos3t-72t2cost) 
+ i(396t cos t - 351 sint - 15 sin 3t + 72t2 sint)] . 



(4.4) 



(4.5) 



(4.6) 



(4.7) 



(4.8) 



Initially prepared coherent state obeys the eigenvalue equation a(0)|a >= a\a >, where a is 
in general complex and the parameter a is given by a = |a| exp(i^) as before. The number 



of photons present in the radiation field (prior to the interaction) is Nq 



\a\ 



For 9 = and 



6 = 7r/2, the corresponding situations are called, in phase (electromagnetic field is in phase with 



the electric field) condition and out of phase condition respectively. Using ( |4-.2[ ) the second order 
variance (hereafter we shall call it as variance) of the X quadrature is calculated with respect to 
the initially prepared coherent state and is given by 



{AXf 



2\a\ 



+ 
+ 



< a\X'^\a > - < a\X\a >^ 

[(|Ei|2 + |E2p(9|«|4 + 18|«|2 + 6) + |^3|'(5|«|^ - 
{2El\a\^a*^ + 2EiEs\a\'^ + EiE*a^ + 3EiE*a*^ 
6E2E3(|a|2a2 ^ ^ 3E2E*a'^ + E^E^a^ + 2EiE^\a\'^a^ 
3EiE*|a|4 + 4EiE*|a|2a*2 ^ ^E^Ela*^ + c.c)] 



(4.9) 



where the terms beyond are neglected. In a similar manner, the variance for other quadrature 
X{t) may be obtained. In absence of the interaction (i.e A = 0), the variance ( |4.9| ) reduces to |. 
The squeezing in absolute sense of X quadrature is obtained if the value of (AX)^ = 0. We note 
that the squeezing of one quadrature automatically prohibits the squeezing of the remaining one. 
The variance (|4.9D is extremely complicated and deserves more investigation for the discussions 
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of squeezing effects of the input coherent light. It is not our purpose to go through such detail 
discussions of the squeezing arising out of the expression (|4.9D. Rather, we would prefer to 
discuss few special cases of physical interests. 

It is true that both the solutions ( 12.25D and ( [2.26| ) are equally good for the purpose of getting 



( AX)^. There are good reasons for using the solution ( [2.25[ ) instead of the solution ( |2.26| ) though 
the previous one contains secular terms. Because of the complex nature of the operators in 
solution ( |2.26| ), it is hard to calculate (AX)^ . In addition to this, the present work takes care 



the solution of the quartic AHO up to the second order in A. Both the solutions ( [2.25D and ( [2.26D 
are identical as long as the second order solution is concerned. The removal of the secular terms 
from the calculated (AX)^, if required, may be done by using the same Tucking-in technique. 

4.1.1 Vacuum field 

In case of vacuum field (i.e |ap = 0), equation ( pr9[ ) reduces to an extremely simple form 

(AX)L = W + 6\E,\' 

= i-^sin2t+ [|fj(201 - 24t2- 208cos2t + 7cos4t-168tsin2t) . 

(4.10) 

If we drop the terms involving present in ( |4.10| ) then the corresponding variance coincides 



exactly with the earlier result [[8q]. Let us now discuss two special cases of particular interests. 
For t = nn, the variance ( |4.10| ) reduces to 



1 QA^ 

(AX)2 =--— ^V. (4.11) 

V /vac 2 

Equation ( |4.1 ID depends on A^ only and we observe that the squeezing of the vacuum field is 
possible for t = nn. Note that a first order calculation can not predict the squeezing effect for 
the same conditions. Hence, a second order solution is not only useful but also essential for the 



present study. For t = f , the variance ( |4.10| ) reduces to 



1 Q \ Q \2 

(AX)L = ^-f + 1^(208 - 3vr^). (4.12) 



For small coupling, the first order term {i.e oc A) in equation ( |4.12D dominates over the second 
order term (i.e oc A^) and hence the squeezing effects are controlled accordingly. Clearly, the 
monotonic increase of the squeezing effects due to the first order term is arrested by the presence 
of the second order one. Thus the second order term is responsible for the reduction of the 
squeezing effect for large coupling. 



Note that the so-called tucking-in technique is also applicable to equation ( |4.10D . Hence, the 
secular terms can be summed up for all orders of A. The corresponding expression for (AX)^^^ 
may be expressed as 

(AX)2^^ = icos2t-f sin2t+g(201 - 208 cos 2t + Tecs 4t) 

+ i(cosf -cos2t(l-^)). ^^-''^ 
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4.1.2 In phase with the electric field 



We assume that the input radiation field is in phase with that of the electric field. In that case 
9 = and a is real. The corresponding variance reduces to the following form 

{AXy = (AX)2„, + ^[-A(l-cos2t + tsin2t) 

+ ^ (11 + 2At^ - 32 cos 2t + 21 cos At + 71t sin 2t + t sin At - 6At^ cos 2t) 

+ [237 + 72t2 - 256 cos 2t + 19 cos At - 36t sin 2t - 36t sin At - 216t^ cos 2t] 

(4.14) 

Equation ( [4-.14D is useful for the purpose of discussing squeezed states generated due to the 



interaction a coherent light with a cubic nonlinear medium. For |ap = 0, corresponding equation 
reduces to the equation ( |4.10D . If we now drop the terms involving A^, the variance ( |4.14D does 



have exact coincidence with the earlier reported results [|6p. Here we are not going to present 
the "Tucked-in form" of the equation ( [4-.14D . This is because, we are interested in the qualitative 
study of (AX)^ up to the second order in A for a weakly coupled system. If we put t = f , the 
corresponding variance reduces to 

(4.15) 



It is obvious that the equation 04-.15D exhibits the monotonic nature of the squeezing as long as 
the first order solution is concerned. In most of the experiments, the percentage of squeezing 
gets saturated after a certain value of the electric field and the squeezing effects are completely 
destroyed if the field strength is very high. Moreover, oscillatory behaviour in the performance 
of squeezing is observed in experiments. These experimental facts are consistent with the present 



work ( |4.15| ). These experimental features also reveals the fact that higher order effects are in- 



evitable in the calculation of squeezing. Thus the present second order calculation is more useful 
to predict the correct behavior of squeezing phenomenon observed due to the interaction of co- 
herent light with dielectric media. 
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Chapter 5 



Photon bunching, antibunching and photon 
statistics 



People tried earnestly to understand the basic nature of the radiation field and as a result of 
this endeavor the concept of photon came into reality [ ]S7| - |S^ . With the due success of the 
photon concept, the immediate questions that come to our mind are: How do the photons come? 
Are they coming together (i.e cluster) or one after another (without cluster)? The answers to 
these questions have started coming with the remarkable experiment of Hanbury-Brown and 
Twiss In their experiment, they obtained the intensity correlation of an incandescent light 
source and they concluded that the photons come together indeed. The phenomenon in which 
the photons come in a cluster is termed as photon bunching. On the other hand, the probability 
of detecting a coincident pair of photons from an antibunched field is less than that from a fully 
coherent light field with a random Poisson distribution of photons. 

The explanations for photon bunching may be given in terms of the wave nature of the light in 
addition to the particle (quantum) nature. However, antibunching of photon can only be explained 
in terms of the particle concept of the radiation field. For this reason, the photon antibunching 
is regarded as a fully quantum mechanical phenomenon without any classical analogue. There 
is a natural tendency that the photons are bunched However, there are several predic- 

tions (models) which go against this tendency and exhibit the photon antibunching. For example, 
various coherent states and the degenerate parametric amplifier [^] were predicted as pos- 
sible sources of getting antibunched light. In addition to this, the resonance fluorescence from 



a two level atom [ |92[1 , the two photon Dick model 0931 1 and the coherent light interacting with 
a two-photon medium [Q] are possible sources for producing antibunched light. The fluores- 
cence spectrum of a single two level atom shows the antibunching and the squeezing as well 
[^2p . There are several examples where the antibunched light are produced in the laboratory [ P5| - 



1O60 . In most of these experiments [ ]97| - |101D , the resonance fluorescence from a small number 



of ions [p7|], atoms or molecules [ PB| - |101[ ] are studied to obtain the antibunching effects. The 
basic physics behind these experiments are clear. The atoms (or molecules) emit radiation and 
go to a ground state from where no subsequent radiation is possible. Hence the emitted photons 
are antibunched. The resonance fluorescence field from many-atom source is not suitable for 
antibunching of photons since the photons emitted are highly uncorrelated. However, a suitable 
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phase matching condition similar to that of four- wave mixing leads to the photon antibunching 
even in the resonance fluorescence of a multi atomic system [ |102| ]. 

With the advent of different type of radiation sources, studies on the quantum statistical 
properties (QSP) of the radiation field are increased considerably. The QSP of the radiation field 
are usually studied with the help of second order correlation function for zero time delay 



107] 



< a)a)aa > 



Equation ( |5.1D can also be expressed in terms of < > and (AA^)^ as 

{ANf- <N > 



9\0) = 1 + 



(5.1) 



(5.2) 



where < > is the average number of photons present in the radiation field and the parameter 
(AA^)2 =< N'^ > — < N >^ is the second order variance in photon number. 

The conditions for bunching and antibunching of photons are given in the table 5.1 



g\0) > 1 


photons are bunched 


g'io) = 1 


photons are coherent 


g\0) < 1 


photons are antibunched 



Table 5.1 



One more useful parameter in the context of QSP of the radiation field is the Mandel Q- 

(5.3) 



parameter [ pO| , |1070 which is defined as 

{ANf- <N > 



Q 



<N>(g\0)-i: 



<N> 

Statistical distribution of photons are related to the Q parameter as 



value of Q 


Photon Number 
Distribution 


Q<0 


sub-Poissonian 


Q = 


Poissonian 


Q>0 


super-Poissonian 



Table 5.2 



Now if we define a quantity 

d = ((AA^)2- < A^ >) , (5.4) 
then the sign of d will essentially determine the quantum statistical properties of the radiation 



field. From the definition of (5.4), Table 5.1 and Table 5.2 we can construct new conditions as 
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Value Ui Ci 


UUllCllillg/ 

fi n ti -Vii 1 n ph i n <T 


r IIUIUII i\ UlllUCl 

F)i 'striViiitinri 





coherent 


Poissonian 


< 


anti-bunched 


sub-Poissonian 


> 


bunched 


super-Poissonian 



Table 5.3 

The purpose of the present chapter is to study the photon number distribution (PND) and 
to investigate the possibilities of getting bunching and antibunching of photons in the interac- 
tion of a coherent light with a third order nonlinear medium of inversion symmetry. We are 
specially interested to keep the off-diagonal terms in the model Hamiltonian and to study the 
consequences of their inclusion in connection with the study of nonclassical effects. The effects 
of input vacuum field on photon bunching and antibunching are also taken care of. 



S.l The photon bunching and photon antibunching 



In this chapter number of photons present in the radiation field prior to the interaction will be 
denoted by nonnegative integer Nq = \a\'^. Therefore, equation ( p.42D may be used to write the 
average number of photons present in the radiation field (initially prepared in coherent state) after 
interaction as 



< a\a^a\a>=< N > = A/'o [l + | {2(3 + 2A/'o) sint sin(t - 2^) 



+ 



(5.5) 



A/'osin2tsin2(t - 2^)}] . 
Thus the number of post interacted photons (i.e < >) is not conserved and is phase sensitive. 



It is not at all surprising since the Hamiltonian ( |1.3lD contains the terms which are nonconserving 
in nature. A considerable amount of discussions about the inclusion and impacts of the photon 
number nonconserving terms in the QAHO model are available in chapter 1 and chapter 3. Now 
the variance (AiV)^ =< a|A^^|Q; > - < a\N\a >^ is given by (see |3.49| ) 



{ANY = A/'o[l + A((3 + 4A/'o)sintsin(t-2^) + A/'osin2tsin2(t-20))], 



Hence we have. 



d 



[(2 + ANo) sin t sin(t - 26) + A^'o sin 2t sin 2(t - 26)] . 



(5.6) 



(5.7) 



The sign of the parameter d controls the phenomena of bunching and antibunching of post- 
interacted photons. However, the oscillatory nature of d ensures the revival and collapse of 
bunching and antibunching effects. The revival of bunching and antibunching was obtained in 
earlier occasion too [|94l]. Now we investigate how the sign of d is controlled by the phases. For 
6 = 0, d is always positive and hence the photons are always bunched. Corresponding PND 
is super-Poissonian. For 6 = f (^), the value of d is negative (positive) if < sin2t < 1 
and is positive (negative) if —1 < sin2t < 0. Thus the photons are antibunched and bunched 
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respectively for the former (later) and later (former) conditions. These statements are further 
corroborated by other values of phases. In the equations ( |53| ) and (|5.6|), the terms beyond linear 
power of A are neglected. This is because of two reasons. Firstly, the solution available for the 



equation ( |1.32D is simple as long as we are interested up to the linear power in A. The second 
reason is that we are interested about the qualitative studies of nonclassical properties of the 
radiation field and hence, the solution linear in A is sufficient. It is tedious, however possible to 
obtain the correction up to the second order in A for d. The second order correlation function 



( ]5^ for zero time delay depends on the photon number Nq, coupling strength (A), free evolution 
time t and on the phase angle 9. Values of t may be controlled by suitable manipulation of the 
length of the interaction region. 

For e = t/2we obtain, < iV >= iVo and (AiV)^ = Nq. Therefore, = and hence the 
input coherent state remains unchanged during its passage through the nonlinear medium. The 
situation may be viewed as the propagation of coherent beam without affecting the propagating 
medium. The particular value of t (= 29) causes a destructive interference in the absorption 
profile and as a result of that we observe a phenomenon which has strong resemblance with the 



self induced transparency (SIT) [|108| 1. Here we can also mention one more important effect 
which is the outcome of the phase of the input coherent light. For 9 = ^ and — 1 < sin2t < 
we obtain the antibunching of the emitted (post-interacted) photons. 

5.2 The effect of vacuum field: 

Now, we give the flavor of analytical results for input vacuum field. The quantized electro- 
magnetic field with no photons are termed as vacuum field. The vacuum field is responsible 
for various interesting physical phenomena. For example, the vacuum field causes Lamb shift, 
Kasimir effects and Spontaneous emission [5^. Note that all these phenomena are purely quan- 



tum electrodynamic in nature and do not have any classical analogue. The effects of vacuum 
field interacting with a cubic nonlinear medium are studied to obtain the squeezing p^] and 



phase fluctuations 01O9| ] of input coherent light. It was found that these nonlinear effects are 



the outcome of inclusion of photon number nonconserving terms. Therefore, it is quite natural 
to explore the effects of vacuum field on the bunching and antibunching of photons. In case of 
input vacuum field {i.e Nq = 0), < N > and {ANY are identically zero and hence the nu- 
merator of the equation ( |5.2D is zero {i.e. d = 0). Hence, g^{0) is indeterminate. However, the 
denominator < faster compare to the numerator. Thus the value of g^{0) ^ 1, if 

the sign of d {i.ed 0+) is positive. Therefore the photons are strongly bunched for vacuum 



field. It is already established that the vacuum field produces the squeezing effect QSq]. Hence 
the retainment of the non conserving energy terms in the model Hamiltonian ( |1.31| ) leads to pho- 
ton bunching and squeezing effects of the vacuum fields. Of course, the two photon coherent 
state with zero photon (i.e squeezed vacuum state) gives rise to the photon bunching along with 



the super-Poissonian photon statistics [90,107]. However the explanation of photon bunching 



for vacuum field were not given in these earlier publications. The situation is interesting, since, 
the photons are bunched although no photons are present in the radiation field! The result does 
not have classical analogue and is a consequence of pure quantum electrodynamic effect. The 
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vacuum field interacts with the medium and produces photons which are bunched. The creation 
of photons through nonlinear interactions is ensured with the incorporation of the photon number 
nonconserving energy terms in the model Hamiltonian ( |1.31D . 

Interestingly, the suitable choice of phases (for example 6 = n /A and < t < 7r/2) leads to 
negative d (i.e d 0^) and hence the photon antibunching. Thus, the antibunched photons may 
be produced through the nonlinear interaction between the medium and the vacuum field. 

For A^o = the Q parameter ( |5.3| ) has indeterminate form (^) . However, in the limit Nq — ^ 0, 



the corresponding value of Q is found to be 

^^■sintsin(t - 29) 
1 + f^sintsin(t - 29) 



Q = . .%A • , • L - (5-8) 



Now for = and t ^ nir, corresponding PND is super-Poissonian. Interestingly, for 6* = n/A, 
the corresponding PND is sub-Poissonian if0<t<7r/2 and super-Poissonian if 7r/2 < t < 
TT respectively. Hence, the sub- and super-Poissonian distribution of the photons are obtained 
periodically. For 9 = n/2, the corresponding PND is sub-Poissonian for all t. 

5.3 Conclusion: 

We obtain the simultaneous appearance of photon bunching (classical phenomenon) and the sub- 
Poissonian photon statistics (nonclassical phenomenon). In other words classical phenomena 
may appear simultaneously with the nonclassical phenomena. These results are consistent with 
the previous observations [ p7| , |l lOQ . The effects of the phase of the input coherent light on g'^{0) 



are pointed out clearly. We obtain the bunching and antibunching of the vacuum field for 9 0. 
It is found that the vacuum field generates photons which may be bunched or antibunched. 
We conclude the present chapter with the following observations. 

1. As a result of the interaction between a single mode of coherent light and a nonlinear 
medium classical and nonclassical phenomena may appear simultaneously. There is no 
reason to believe that the appearance of one of the classical (nonclassical) phenomena 
warrants the presence of the others. 

2. We report the sub-Poissonian PND and the antibunching of photons for input coherent 
light and vacuum field. 

3. The vacuum field interacting with the nonlinear medium produces photons through non- 
linear interaction. 
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Chapter 6 

Application of the m-th anharmonic 
oscillator: Interaction of coherent light 
with an (m — l)-th order nonlinear medium 



The concept of phase plays a very crucial role in the understanding of basic physics. In fact, phase 
is responsible for all the interference phenomena we observe in classical and quantum physics. 
So people have studied the properties of the quantum phase from different contexts since the 
beginning of quantum mechanics [ |111| ]. But the interest increased considerably in the recent past 
when Berry published his famous paper []67|]. He showed that the phase acquired by the quantum 
system during a cyclic evolution under the action of an adiabatically varying Hamiltonian is a 
sum of two parts. The first part is dynamical in nature and the second is geometric. Later on 
Aharonov and Anandan generalized the idea of the Berry phase and they defined a geometric 
phase factor for any cyclic evolution of quantum system [|68|]. The existence of geometric phase 
is found in many areas of physics. The examples are from the quantum Hall effect to the Jahn- 
Teller effect and from the spin orbit interaction to quantum computation P Jl 1 1|JTT2| ]. 

In the present work, interaction of electromagnetic fields with matter is modeled by AHO. 
These facts provoked us to study the possibilities of observing nonadiabatic geometric phase or 
Aharonov-Anandan phase [ [6^J1 13| ] for anharmonic oscillators in general. In the present chap- 
ter we will find an analytic expression in closed form for the Aharonov-Anandan phase for a 
generalized anharmonic oscillator. The generalized expression is then used to study a particular 
case of physical interest in which an intense laser beam interacts with a third order nonlinear 
nonabsorbing medium having inversion symmetry. In this chapter results of chapters 4 and 5 are 
extended to the case of (m — l)-th order nonlinear medium in general to study the possibilities 
of generating nonclassical states in an (m — l)-th order nonlinear medium. 



6.1 Aharonov-Anandan phase 

In quantum mechanics of any system, one has a complex Hilbert space H. whose nonzero vectors 
represent states of the system. A vector > mH and any complex multiple of it (i.e X\il} >) 
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represent the same state in the Hilbert space H. This arbitrariness in the representation of states 
can be reduced by imposing the normalization condition 



<V|^>=1- (6.1) 

There still remains an arbitrariness of a phase factor because a normalized vector \ip > and the 
vectors a\4' > represents the same state in the Hilbert space provided that the modulus of a being 
unity (i.e \a\ = 1). The collection of all vectors > with |^ >, a fixed normalized vector 
and a taking all possible complex values of modulus one is called a unit ray and is denoted by 
IV' > [ |1 14| ]. More general rays are consist of collection of all vectors of the form X\tp >. So rays 
represent vectors without any arbitrariness of phase. 

To understand the geometric phase let us start from a simple situation in which a state vector 
\ip > evolves cyclically in the Hilbert space H and after a complete period T returns to the same 
state vector with a different phase which may be written as a I?/) >= exp(i0) {i/j >. Therefore, 

>=exp(z0)|^(O) > (6.2) 

where (p is the total phase which has two parts. The first part is the dynamical phase and the other 
part is the nonadiabatic geometric phase or the Aharonov Anandan phase (AA phase). We are 
interested in the later type of phase which can be found out by subtracting the dynamical phase 
part from the total phase. 

The cyclic evolution of the state vector may be described by a curve C inH because \ip > and 
a\ip > are two different points in H. Since there are infinitely many possible values of a so there 
are infinitely many curves C inH which describes the same cyclic evolution. In the projective 
Hilbert space of rays P or in the ray representative space all these points are represented by the 
same ray (i.e the same point). So we will have a single closed curve C in the ray representative 
space P corresponding to the infinitely many possible curves in the usual Hilbert space. As there 
are no arbitrariness of phase in P so in the interval [0, T] we must have 

>= \ij{0) > . (6.3) 

Aharonov and Anandan exploited this property of the single valuedness of the unit ray in the ray 
representative space P to establish the presence of a geometric phase for all cyclic evolutions. 
They started their work by supposing that the normalized state >E'H evolves according to 
the Schrodinger equation 

H{t)\m>=^f^[j^\m>y (6.4) 

We can always define a state vector \^{t) > in P which is equivalent to \ilj{t) > as 

m)>=exp{-zfmm> ■ (6.5) 



Now substituting ( p3 and ( |63| ) in ( |63| ) we obtain 

/(T) - /(O) = (6.6) 
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If we assume that the total wave function rotates by 2tt radian for this T which results in a 
cyclic motion of every state vector of the Hilbert space H, then we obtain another condition as 



/(T) - /(O) = 2n. (6.7) 

Mere satisfaction of condition ( |6.7D will ensure the single valuedness of the wave function in the 
usual Hilbert space H. Again from ( |6.4| ) and (^3|) we have 

H\m >= (J-^exp{zf{t)m) >) (6.8) 

or, 

H\ij{t) >= -h^^exp{tf{t)m) > +ihexp{zf{t)j^m) > . (6.9) 

Therefore, 

< m\H\m >= -^f < ^(t)|exp(z/(t)|^(t) > +th < ^(t)|exp(^/(t)^|^(t) > 

(6.10) 

or, 

- f < m\m >=l< m\H\m >-< mi^j^im > ■ (6.11) 

Since ijj{t) is normalized so we have 

- f = ^ < ^(^)\H\m >-< > . (6.12) 

at n at 

Now integrating equation ( |6.12p from to T and using ( |^^ we have 

^ < MIHIM) >dt+ r < Mli^M) > dt = -f + (3 (6.13) 

JO at 



7 = -^ r <m\H\m>dt (6.14) 
n Jo 



where 

-T 

h Jo 

is proportional to the time integral of the expectation value of the Hamiltonian, it is dynamical in 
origin and hence called dynamical phase. Now if we subtract the dynamical phase part 7 from 
the total phase then we get the geometric phase or the Aharonov Anandan phase 

(3= r <i,{t)\i^Mt)>dt. (6.15) 
JO at 

It is clear that the same > can be chosen for infinitely many possible curves C by ap- 

propriate choice of f{t). This phase is geometrical in the sense that it does not depend on the 
Hamiltonian responsible for the evolution. Moreover it is independent of the parameterization 
and redefinition of phase oi\il){t) >. Hence, P depends only on the evolution of the shadow of 
\^{t) > in projective Hilbert space P and is a pure geometrical entity. 
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6.1.1 Anharmonic oscillator: 



The Hamiltonian of a generalized anhamionic oscillator can be written as 



H = Ho + Hi 

= a^a + X'F{a^ , a) 



(6.16) 



Now the total eigenstate of the unperturbed Hamiltonian Hq is 

\^{0) >=Y,Cn\n> (6.17) 

n 

where the coefficients Cn's depend on the initial conditions. The total wave function of the 



Hamiltonian H satisfying the Schrodinger equation ( |6.4[ ) is 

lipit) >= ^ Cn exp (y—iHt) \n >= ^ C„ exp ( — ia^at — i}>! F{a\ a)t] \n > . 



(6.18) 



For the calculation of the geometric phase, we have to chose the single- valued state >. 
This can be easily done by choosing a f{t) which simultaneously satisfies equation ( |01 ) and 



(p7l\). The simplest choice is f{t) = X't and with that choice we have 



m)> = expi-tfmm> 



J2n exp (^—ia^at — iX't (l + F{a'', a 



\n > 



(6.19) 



and 



/(T) = X'T = 271. (6.20) 

Now from equations ( |6.15| , |6.19| and |6.20D a compact expression for the Aharonov Anandan 
phase of a generalized anharmonic oscillator is obtained as 

p = <m\ti\m>dt 

= lo Efc En< k\ClCn exp [ia^at + iX't {l + F{a\ 
X a) a + A' (l + F{a\ a)) exp [-a) at - iX't (l + F{ 
= lo En < k\ClCn [a^a + A' (l + F(at, a) ^ 



\n > dt 



n > dt 

n > 



(6.21) 



= EfcEn< A;|QC„[27r(l + ^ + F(at,, 

= 2tt + En \Cn\^ (^) + 2vr Ek En < (F(at, a)) \n > . 

The quantum nature of the geometric phase is manifested through the discrete sum appearing 
in the expression ( |6.21D and the state dependence of the phase appears through the coefficients 
C^Cn- Now let us think of some special situations of physical interest. 



6.1.1.1 Anharmonic part is a polynomial of the number operator: 

When the anharmonic part can be written as a polynomial P of number operator a^a, i.e when 
F{a\ a) = P{a)a) then from ( |6.21| ) we obtain 



/? = E|Cn| 



27r 1 + 



n 
J' 



Pin 



(6.22) 



79 



One of the simplest possible polynomial of number operator is 



2 

a' a 



a'' a. 



This kind of interaction is very common in Physics and Joshi et al [0] have shown that a geo- 
metric phase exists for such a Hamiltonian. Actually they considered the propagation of elec- 
tromagnetic radiation in a dispersive medium and under certain approximation they obtained the 
approximated Hamiltonian. Now from ( |6.22| ) we obtain the geometric phase for the approxi- 



mated Hamiltonian of Joshi et al (1.33) as 



n 



27r ( 1 + - + (n^ 



n 



(6.23) 



This expression for geometric phase does not coincide exactly with that of Joshi et al [^. The 
term En |Cn| 



27r^ 



is absent in their calculation. A conceptual oversight is responsible for 
that. In their paper they have worked in the interaction picture and found out the correct and 
exact wave function (in the interaction picture) but they have not considered the time evolution 



of the differential operator present in equation 



).15D . The time evolution of the operator 



dt 



in the interaction picture is exp{iHot)-^ exp{—iHot) and this is the operator one should use in 



equation ( |6.15D while working in the interaction picture. Use of the corrected operator will yield 
the same result. 



6.1.1.2 Interaction with an (m-l)-th order nonlinear medium 

Now let us think of a physical situation in which an intense electromagnetic field interacts with 
an (m — l)-th order nonlinear nonabsorbing medium. If the symmetry of the medium is chosen 
in such a way that the only nonlinear interaction in the medium appears due to the presence of 
(m — l)-th order susceptibility then the interaction Hamiltonian ( |1.29| ) of the system may be 
written in an alternative form as 

Tj _ ^ _i_ ^ _\ A_ m 

~ 2 2 m-^ (f. 24-') 

= ata + A'(at + a)™ ^ 



where A' = ^^^m . The above Hamiltonian ( |6.24| ) represent a generalized anharmonic oscillator 



for which F(a^, a) = (a^ + a)"^ and the geometric phase is 

(3 = 27r + E„|C„pfe)+27rEfeE„Cfe*C„<A;|(at + a)-|n> 

= 27r + En ICnP (^) + 2vr En ErU T.p=o' hr "^"^^C^ (6.25) 

^* [n!(ra-m+2r+2p)!]^ 



where 



2(5)(i)! 



(6.26) 
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Normal ordering theorems 1 and 2 are used here to derive (|6.25|). Now we have a closed form 



analytic expression ( |6.25D for the geometric phase (3 for a generalized anharmonic oscillator 
which depends on the photon statistics of the input radiation field. So an observation of the 
geometric phase may help us to conclude the nature of the input electromagnetic field. We use 
this generalized expression ( |6.25| ) to study a very special case. 

6.1.1.3 An intense laser beam interacts with a third order nonlinear medium 

Let us now consider the case where an intense electromagnetic field having Poissonian statistics 
(C„ = exp(— ■^)^) interacts with a third order nonlinear nonabsorbing medium. If the sym- 
metries of the medium is chosen in such a way that the only nonlinear interaction in the medium 
appears due to the presence of third order susceptibility, then the Hamiltonian of the system is 

^ = y + y + = ^0 + A'(at + af (6.27) 



where A is the coupling constant and A' = ^. With the help of equations ( |6.25D the geometric 
phase [3 for this physical system can be written as 

a — 9^_i_V ir* |2 /^27m\ I 9 v^4-2r^ iri i-2r [n!(n-4+2r+2p)!]^ 

= 87r + 27rexp(-|«n{E„M!:i [(Gn^ + (6 + i) n) 
= 4|a|2 {2n + 3) cos(2^) + 2|a|2 cos(4^)]} 

= Stt + 27r [(6|a|^ + 12|a|2 + ^) + 4|a|2 {2\a\^ + 3) cos(20) + 2|a|2 cos(4^)" . 

(6.28) 

Here we observe that the geometric phase depends on the phase of the initial laser beam. The 
presence of the off-diagonal terms in the interaction Hamiltonian is manifested through the pres- 
ence of 9 dependent terms in the expression for the geometric phase p. 

6.1.2 Remarks on the results: 

The geometric phase for the physical systems considered by us can be experimentally observed 
with the help of interferometric experiments. To be more precise we can think of an optical 



analogue of the Suter-MuUer- Pines experiment Ql ISQ . A gedanken experiment for this purpose 
is also proposed by Joshi et al. [0]. 

We have seen that the geometric phase arising due to the interaction of intense electromag- 
netic field with a nonlinear medium is experimentally observable and the C„'s depend on the 
photon statistics of the input field (see Table 6.1). So an observation of the geometric phase can 
give us the photon statistics of the input electromagnetic field. Earlier we have shown that the 
photon statistics in a third order nonlinear nonabsorbing medium of inversion symmetry changes 
with the interaction time and in this chapter we will show that this fact is true for {m — 1) — th 
order nonlinear medium also. Using the output of such an interaction as the input of a geometric 
phase experiment our predictions on photon statistics can be verified. 
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Photon Statistics 




Sub-Poissonian 


jj^ri^-pr- {o<p<i) 


Poissonian 


exp(-|a|^)^ 


Super-Poissonian 





Table 6.1 



The general Hamiltonian of the form ( |6.24D and its special case ( |6.27D are very common in 
quantum optics and we often use rotating wave approximation (RWA) to deal with these Hamil- 
tonians. Under RWA calculations off-diagonal terms are neglected. But in the present case the 
presence of the off-diagonal terms in the interaction Hamiltonian is manifested through the pres- 
ence of 9 dependent terms in the expression for the geometric phase p. From these facts we can 
conclude two things firstly, it is not right to use RWA for the calculation of the geometric phase 
and, secondly, geometric phase can be tuned because it depends strongly on the phase of the 
input field which can be tuned. 



6.2 Higher harmonic generation 

If an input of frequency cu generates a output frequency nu due to the nonlinear frequency mixing 
in a nonlinear medium, then we say that the n — th harmonic generation is taken place in the 
medium. It is well known that the n — th harmonic generation is observed in an n — th order 
nonlinear medium. Now, a close look at the term, exp {iVLt{2p — m + 2r — 1)) in ( |2.8(J| ), will 
show that the expansion of the summation present in ( |2.80D will give us terms with frequency 
3f2, 5fi, 7f2, (m — 1)VL when m is even and 2fi, AVt, 6fi, (m — 1)VL when m is odd. Thus 



we can conclude from the above discussion and ( |2.80D that the frequency mixing occurred in 
such a way that all the odd / even higher harmonics up to the order of the nonlinearity of the 
medium will be present in the output. 



6.3 Bunching, antibunching and statistical distribution of the 
photons 

In chapter 5 we have seen that the quantum statistical properties of radiation field can be studied 
with the knowledge of sign of d = ((AA^)^— < >). To study the QSP of the output of an 
interaction of coherent light with an (m— l)-th order nonlinear medium we can construct a closed 
form expression for d by using ( |2.78[ ) and ( |5.4D as 



d 



4A 



X sm 



■ ^ , ^ V, tor 

(2p-m+2r) 
2 



lal™ ^'"sin 



(2p-m+2r) 
2 



(t - 29) 



(6.29) 

where we have taken the expectation value with respect to the initial coherent state \a >= 
X^J^o^^P (^) T^il^ ^ ^^'•^ ^ ~ |q^| exp(z6'), having Poissonian statistics. From this ex- 
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pression we can see that for m = 4, i.e for a third order nonlinear medium having inversion 
symmetry, we have 



d 



3A|a| 



2 (2|ap + l) sin(t - 26) sin(t) + | sin (2(t — 29)) sin(2t) . (6.30) 



This is in exact accordance with our earlier result. From (|6.29|) we have following observations, 

i) When t = 26 then = 0. Therefore, input coherent state remains coherent in the output. 

ii) When 6* = or 6* = nvr, i.e. input is real then d is a sum of square terms only. So d is 
always positive and we have bunched photons having super-Poissonian statistics in output. 

iii) For other values of phase 6 of input radiation field, value of d oscillates from positive 
to negative, so we can observe bunched, antibunched or coherent output depending upon the 
interaction time t. 



6.4 Squeezing 

Here we will study the possibilities of electrically squeezed field. Possibilities of magnetically 
squeezed field may also be studied in the similar process. Now the uncertainty in electric field is 



[AXf 



1 



2\t v^T 



2^m 



E2 + m,r~< 
r=0 ''Zr <-^: 



2^ \-y m — 2r 



' m—2r \ I m—2r 



2 'y 2 

|m — 2r-2 



m-2r-2 



sin {2{e-t)) 



- Z-*" ^2r (-.2r Z.p_^m^ ' ^P(2p-m+2r) \ 2 



2^m 



X 



{p — 1) sin 



{m-2r-2p+2) 



{29 -t)-t) + {m-2r- 2p) sin ( ("^-2j-2p) (2g _ t)) 



(6.31) 

Here we can note that the squeezed state may be generated due to the interaction of an intense 
laser beam with the (m — l)-th order nonlinear medium for particular values of 9 and t. For 
example, for = and t < ^ we will get electrically squeezed field in the output. To compare 
the above expression with the existing results we put m = 4 and 9 = 0, which corresponds to the 
case in which an intense laser beam having zero input phase interacts with a third order nonlinear 
non absorbing medium of inversion symmetry and obtain, 

(AX)2 = ^ - Y si^'(^) - + 2 sm\t)) . (6.32) 



This is in exact accordance with our earlier result []3S[]. From ( |6.32p we can see that the squeezing 
in electric field quadrature will be observed for t < | for any arbitrary allowed values of other 
parameters. So a nonclassical phenomenon (squeezing) is observed for m = 4 and = but 
for the same conditions photons of the output radiation field are always bunched. Thus we can 
conclude, although the nonclassical phenomenon of antibunching and squeezing usually appears 
together but it is not essential that the two nonclassical phenomena have to appear together. 
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Chapter 7 

Summary and concluding remarks 



In the present work coherent Ught interacting with a nonUnear medium of inversion symmetry is 
modeled by a general quantum anharmonic oscillator. The quartic-anharmonic-oscillator model 
emerges if only the lowest order of nonlinearity (i.e third order nonlinear susceptibility) is to be 
present in the medium. The quartic-oscillator model has considerable importance in the study 
of nonlinear and quantum optical effects that arise in a nonlinear medium of inversion symme- 
try. For example, silica crystals constitute an inversion symmetric third order nonlinear medium 
and these crystals are used to construct the optical fibers. In the optical communication electro- 
magnetic beam passes through the optical fiber and the interaction of the electromagnetic field 
(single mode) with the fiber can be described by the one dimensional quartic anharmonic oscil- 
lator Hamiltonian. Depending on the nature of nonlinearity in a physical problem the treatment 
of higher anharmonic oscillators also assumes significance. But anharmonic oscillator models 
are not exactly solvable in a closed analytic form. On the other hand, we need operator solu- 
tions of the equations of motion corresponding to these models in order to study the quantum 
fluctuations of coherent light in nonlinear media. So we have two alternatives, either we can use 
an approximate Hamiltonian which is exactly solvable or we can use an approximate operator 
solution. Approximate operator solutions of anharmonic oscillator problems (solutions in the 
Heisenberg approach) were not available even in the recent past presumably because the existing 
methods tend to introduce inordinate mathematical complications in a detailed study. Due to 
the unavailability of the operator solution people would use RW approximated Hamiltonians to 
study the quantum fluctuations of coherent light in nonlinear media. But the situation is now im- 
proved and many proposals to obtain approximate operator solutions of anharmonic oscillators 
have appeared. Some of these solutions are obtained as a part of the present work. For example, 
we have constructed a second order operator solution for quartic oscillator and have generalized 
the first order operator solutions available for the quartic oscillator to the m-th anharmonic oscil- 
lator. From the generalized solutions we have observed that an apparent discrepancy is present 
between the solutions obtained by different techniques. Then the question arise: Which solution 
should be used for physical applications? Therefore, we have compared different solutions and 
have concluded that all the correct solutions are equivalent and the apparent discrepancy what- 
soever is due to the use of different ordering of the operators. These solutions are then exploited 
to study the possibilities of observing different nonlinear optical phenomena in a nonlinear di- 
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electric medium. To be precise, we have studied quantum phase fluctuations of coherent light 
in third order inversion symmetric nonlinear medium in chapter 3. Fluctuations in phase space 
quadrature for the same system is studied in chapter 4 and the possibility of generating squeezed 
state is reported. In chapter 5 fluctuations in photon number is studied and the nonclassical phe- 
nomenon of antibunching is predicted. In chapter 6 we have generalized the results obtained 
for third order nonlinear medium and have studied the interaction of an intense laser beam with 
an (m — l)-th order nonlinear medium in general. Aharonov Anandan nonadiabatic geometric 
phase is also discussed in the context of (m — l)-th order nonlinear medium in general. 

We have observed different interesting facts regarding quantum fluctuations of coherent light 
interacting with nonlinear media and now we can conclude the present work with the following 
observations. 



1 . In a third order nonlinear medium reduction and enhancement of phase fluctuation param- 
eters (U {9, t, lap), Q {6, t, lap) and S {6, t, |a;p)) are possible for suitable choice of free 
evolution time t. This result is in sharp contrast with the earlier results [^^. 

2. Electrically squeezed electromagnetic field can be generated due to the interaction of an 
intense beam of coherent light with an (m — l)-th order nonlinear medium. 

3. Antibunching of photons can be observed due to the interaction of an intense beam of 
coherent light with an (m — l)-th order nonlinear medium in general. 

4. Simultaneous appearance of photon bunching (classical phenomenon) and squeezing (non- 
classical phenomenon) is observed. This establishes the fact that there is no reason to be- 
lieve that the appearance of one of the nonclassical phenomena warrants the presence of 
the others. 

5. The vacuum field interacting with the nonlinear medium produces photons through non- 
linear interaction. We report sub-Poissonian PND and antibunching of these photons. 

6. A nonvanishing noadiabatic geometric phase appears due to the interaction coherent light 
with an (m — l)-th order nonlinear medium. The geometric phase for the physical sys- 
tems considered by us can be experimentally observed with the help of interferometric 
experiments. 

7. Under RWA calculations off-diagonal terms are neglected. But in the present case the 
presence of the off-diagonal terms in the interaction Hamiltonian is manifested through 
the presence of 9 dependent terms in the expressions for the P, (AX)^ etcetera. Therefore 
special care should be taken before using RWA for the calculations related to the matter- 
field interaction. 



8. Percentage of squeezing and geometric phase can be tuned because they strongly depend 
on the phase of the input coherent field which can be tuned. 
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7.1 Limitations and scope for future works 



In most of the studies envisaged in the present work we have used either first order or second 
order solutions. But higher order terms might have important effects in comparatively strongly 
coupled systems. Higher order operator solutions are now appearing in the literature such that it 
is now technically possible to study the strongly coupled systems using the methods used in the 
present work. 

In some particular cases we have used secular solutions to study the quantum fluctuations. 
This is valid for short interaction time. In order to have corresponding expressions valid for all 
times we have to remove the secular terms from the expressions. Tucking in technique may be 
used for the purpose. 

We started working on the quantum anharmonic oscillator problem in 1998. At that time only 
two techniques were available to provide operator solution of the quantum quartic anharmonic 
oscillator. The situation has considerably improved during the last four years. Now we have sixth 
order solution of quartic oscillator [ |1 16| ]P]- So we can now address the question of convergence of 
the operator series. We can also extend our works on one dimensional anharmonic oscillators to 
two or higher dimensions. We also expect that the present work will serve as the back bone of the 
nonlinear dynamical studies (in Heisenberg approach) of more complex systems. For example, 
we can study a system of coupled oscillators, toda lattices etcetera. This kind of study may 
provide many new information in near future. 

The other thing is that almost all the effects discussed in the present thesis are experimentally 
realizable and all of them have potential applications in optical communications and other areas 
of physics which are related to public life. We finish the present work with an expectation that 
it will be useful for the future development of theoretical and experimental studies of nonlinear 
dynamical systems. 



' In fact we can construct solution of anharmonic oscillators (for any particular m) up to arbitrary order The 
procedure is reported in our work [ |40| ] but it is not included in the present thesis. 
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